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Abstract 


Liquid-vapor (bubble) interface disturbances caused by various types of 
accelerations, including centrifugal, lateral and axial impulses , gravity 
gradient and g-j itter accelerations associated with spinning and slew motion in 
microgravity, have been reviewed. Understanding of bubble deformations and its 
fluctuations are important in the development of spacecraft orbital and attitude 
control techniques to secure its normal operation. This review discusses bubble 
deformations and oscillations driven by various forces in microgravity 
environment. The corresponding bubble mass center fluctuations and slosh 
reaction forces and torques due to bubble deformations have also been reviewed. 
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I. Introduction 


For the purpose to carry out scientific experiments, some experimental 
spacecraft use cryogenic cooling for observation instrumentation and telescope, 
superconducting sensors for gyro read-out and maintain very low temperature near 
absolute zero for mechanical stability. The approaches to both cooling and 
control involve the use of superfluid liquid helium II. In this work, sloshing 
dynamics associated with spinning and/or slew motions are reviewed and 
investigated. The potential problems for cryogenic liquid in dewar container 
could be due to asymmetry in the static liquid helium distribution and to 
perturbations in the liquid-vapor interface (bubble) caused by slosh wave 
excitation driven by pointing control, machinery vibration, etc. 

With liquid helium II, at a temperature below a A-point (2.17 K) , there is 
negligibly small temperature gradients. In the absence of temperature gradient 
along the surface which drive Marangoni convection, the equilibrium shape of the 
free surface is governed by a balance of capillary, centrifugal and gravitational 
forces. Determination of liquid-vapor interface (bubble) profiles based on 
computational experiments can uncover details of the flow which can not be easily 
visualized or measured experimentally in a microgravity environment. 

The instability of the liquid-vapor interface surface can be induced by the 
presence of longitudinal and lateral accelerations. Slosh waves are, thus, 
excited which produces high and low frequency oscillations in the liquid systems. 
The sources of the residual accelerations range from the effects of the Earth's 
gravity gradient 1-3 and g-jitter 4 accelerations which include, atmospheric drag 
on the spacecraft, vibration of compressor, spacecraft attitude motions arising 
from machinery vibrations, thruster firings, spacecraft slew motion, pointing 
control of spacecraft, crew motion, etc. A recent study 5 suggests that the high 
frequency accelerations may be unimportant in comparison to the residual motions 
caused by low frequency accelerations. 
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Time-dependent dynamical behavior of partially-filled rotating fluids in 
reduced gravity environments was simulated by numerically solving the Navier 
Stokes equations subject to the initial and boundary conditions 6-20 . At this 
interface between the liquid and the vapor fluids, both the kinematic surface 
boundary condition, and the interface stress conditions for components tangential 
and normal to the interface, were applied 10-15 . The initial conditions were 
adopted from the steady-state formulations developed by Hung et al 7-9 . Some of 
the steady— state formulations of interface shapes were compared with the 
available experiments carried out by Leslie 21 in a free-falling aircraft (KC- 
135). The experiments carried out by Mason et al 22 showed that the classical 
fluid mechanics theory is applicable for cryogenic liquid helium with 
sufficiently large velocities, exceeding critical velocity, and also in the 
larger containers 23-25 . As to the dynamical behavior of cryogenic helium bubble 
deformations in microgravity, there is only very limited numbers of poor quality 
record for the scientific experimental observation available in a very short time 
interval and in reduced but not microgravity environment 26 , . It is anticipated 
that microgravity experiments will be carried in the near future shuttle flights 
by NASA scientists. In the computer code validation, results of model 
computation for cryogenic liquid draining and shut— off in reduced gravity with 
geyser excitation have been compared with experimental observation with excellent 
agreement 27 . However, the results of mathematical modelling is still urgently 
important to provide the key data for the design of the scientific spacecraft , 
in particular, for their design of guidance and attitude control systems. 

During the spacecraft orbit around the Earth, the direction of azimuth 
angle of Earth toward the location of the spacecraft mass center varies from 0° 
along the spin axis of spacecraft to a full circle of 360° which requires the 
adaptation of three dimensional calculation. 

As the spacecraft moves along the orbit, any fluid capable of motion 
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relative to the spacecraft is subject to the acceleration that arises from the 
gravity gradients of the Earth which is essentially tidal in nature. Detailed 
description of the characteristics of gravity gradient will be given in 
Section(II-D) . The interaction between the mass of fluid element and the 
spacecraft mass due to gravity gradient accelerations 1 " 3 are capable to excite 
slosh waves and induce the fluctuations of slosh reaction forces and its moment 
acting on the dewar of spacecraft fluid systems. Slosh wave excitations in a 
spinning container of a linearized fluid have been carried out by Bauer 28 * 29 . 

At temperatures close to absolute zero, quantum effects begin to be of 
importance in the properties of fluids. At a temperature of 2.17 K, liquid 
helium has a A-point (a second-order phase transition) ; at temperatures below 
this point, liquid helium (helium II) has a number of remarkable properties, the 
most important of which is superfluidity. This is the property of being able to 
flow without viscosity in narrow capillaries or gaps. At temperatures other than 
zero, helium II behaves as if it were a mixture of two different liquids. One 
of these is a superfluid and moves with zero viscosity along a solid surface. 
The other is a normal viscous fluid. The two motions occur without any transfer 
of momentum from one to another for velocities below a critical velocity 23-25 . 
For the components of normal and superfluid velocities above a critical velocity, 
the two fluids are coupled 23-25 . 

Experiment made by Andronikashvili 23 * 24 for rotating helium II shows that 
it is necessary to exceed a critical velocity for the interaction between the 
normal and superfluid components to establish entire bucket in rotation 23 * 24 . 
This is due to the fact that the high fluid velocity can produce great enough 
vortex lines to snarl in a complex tangle to assure and interaction between the 
normal and superfluid components 23-25 ' 28 * 29 . For the rotating dewar with outer 
diameter of 1.56 m and Inner diameter of 0.276 m, the critical velocities to 
assume the interaction between the normal and super— fluid components are 6.4xl0 -7 
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and 3.6xlCT 6 m/s, respectively 23,2 * . With rotating speed of 0.1 rpm, the linear 
velocities along the outer and inner walls of rotating dewar are 8.17xlCT 3 and 
1.45xl0“ 3 m/s, respectively, which are at least several hundred times greater 
than that of the corresponding critical velocities to assure the interaction 
between the normal and superfluid components of helium II. These vortex lines 
snarling with a complex tangle between the normal and superfluid components 
warrant the adoption of Newtonian fluid model 23-25 . Based on this illustration, the 
problem under consideration have the special features to warrant an adoption of 
viscous Newtonian fluids formulation in this study. 

Density concentration of superfluid is a function of temperature, which 
is also true for the surface tension and viscous coefficient for helium II 30,31 . 
In general two fluid model shall be adopted to compare the dynamical behavior of 
helium II. For the cases with fluid velocities exceed a critical velocity, 
single fluid model can be employed. In this study, the theory of viscous 
Newtonian fluids is employed with modification of transport coefficients adjusted 
by normal and superfluid density concentration which is a function of 
temperature. Most of the research works pertaining to the cryogenic helium 
sloshing dynamics in microgravity are carried out by NASA Marshall Space Flight 
Center (NASA/MSCF) , while some of the cryogenics hydrogen sloshing dynamiics in 
reduced gravity are operated by NASA Lewis Research Center 26 . Majority of the 
present works of cryogenics helium II are based on the works accompanied by a 
group of scientists, affiliated with the NASA/MSCF in the past 8 years. In 
addition to the application of helium II sloshing dynamics with velocities 
exceeding critical velocity in microgravity, present works are also fully 
applicable to sloshing dynamics with regular Newtonian fluids. 

For the purpose of examining the sloshing dynamics affected liquid-vapor 
interface in a partially liquid-filled dewar subject to various conditions of 
spacecraft orbiting around the Earth, examples are given to illustrate the 
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dynamics of bubble deformations under the following conditions: (1) Bubble spin- 

up from rest with and without a completion of wrapping around the well of 
rotating dewar, (2) Bubble spin-down from steady state with wrapping around 
central column of rotating dewar, (3) Rotating bubble subject to lateral and 
axial impulses, (4) Bubble subject to gravity gradient or g-jitter accelerations 
associated with slew motion, and (5) Rotating bubble subject to various 
magnitudes of combined gravity gradient and g-jitter accelerations. 

Dynamics of Bubble deformation introduce additional problems of spacecraft 
orbital mechanics in guidance and attitude controls under microgravity 
environment. These problems include: (a) bubble mass center fluctuations, and 

(b) fluctuations of slosh reaction forces and torques exerted on the liquid 
container of spacecraft fluid systems due to bubble deformations. In this review 
paper, it is our intention to discuss (a) the cross-linkage between bubble 
deformation driven by various orbital accelerations and their impacts on bubble 
mass center fluctuations, and (b) the corresponding fluctuations of slosh 
reaction forces and torques acting on the fluid container due to bubble 
deformations. Better understanding of these linkages is a key to develop the 
reliable spacecraft guidance and attitude control techniques. 

Present research has resulted in numerous refereed journal article 
publications. List of journal article publications sponsored by present NASA 
Grant Number: NAG8-938 are as follows: 

(1) Hung, R. J., and Shyu, K. L. , Constant Reverse Thrust Activated Reorientation of Liquid Hydrogen with 
Geyser Initiation, Journal of Spacecraft and Rockets . 29 . 279-285, 1992. 

(2) Hung, R. J. , and Shyu, K. L. , Fluid Resettlement and Slosh Wave Excitement During the Spacecraft 
Reorientation in Cryogenic Propellant System, Proceedings of National Science Council (A), 1J5, 131-148, 
1992. 

(3) Hung, R. J., Lee, C. C., and Leslie, F. W. , Spacecraft Dynamical Distribution of Fluid Stresses 
Activated by Gravity Jitters Induced Slosh Waves, Journal of Guidance. Control and Dynamics . 15 . 
817-824, 1992. 

(4) Hung, R. J . , Shyu, K. L., Lee, C. C . , Liquid Hydrogen Slosh Wave Excited by Constant Reverse Gravity 
Acceleration of Geyser Initiation, Journal of Spacecraft and Rockets , 29, 523-528, 1992. 

(5) Hung, R. J., Shyu, K. L. , and Lee, C. C., Medium Frequency Impulsive Thrust Excited Slosh Waves 
During Propellant Reorientation with Geyser, Journal of Propulsion and Power . 5* 778-785, 1992. 

(6) Hung, R. J. , and Shyu, K. L.. Excitation of Slosh Waves Associated with Low Frequency Impulsive 
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Reverse Gravity Acceleration of Geyser Initiation, Acta Astronautics . 26, 425-433, 1992. 

(7) Hung, R. J., and Lee, C. C. , Characteristics and Behaviors of Gravity Probe-B Spacecraft Propulsion 

System, Proceedings of National Science Council (A), lj>, 239-252, 1992. 

<8 ) Hung, R. J. f and Shyu, K. L., Suction Dip, Liquid Residual and Liquid Propellant Shut-Off in 

Spacecraft Propulsion System, Proceedings of National Science Council (£) , ^6, 347-364, 1992. 

(9) Hung, R. J., Lee, C. C. , and Leslie, F. W. , Similarity Rules in Gravity Jitter-Related Spacecraft 
Liquid Propellant Slosh Waves Excitation, Journal of Fluids and Structures . £, 493-522, 1992. 

(10) Hung, R. J., and Shyu, K. L. , Medium Frequency Impulsive Thrust Activated Liquid Hydrogen 
Reorientation with Geyser, Journal of Propulsion and Power . £, 987-994, 1992. 

(11) Hung, R. J. , Lee, C. C. , and Leslie, F. W,, Effect of the Baffle on the Spacecraft Fluid Propellant 
Viscous Stress and Moment Fluctuations, Transactions of the Japan Society for Aeronautical and 
Space Sciences . 35, 187-207, 1993. 

(12) Hung, R. J., and Shyu, K. L., Suction Dip, Liquid Residual and Slosh Wave Excitation During Liquid 

Draining in Microgravity, Advances in Space Research . 13(7) . 147-154, 1993. 

(13) Hung, R. J. , and Shyu, K. L., Liquid Settlement, Resettlement, Slosh Wave Excitation and Geyser 

Motion During Reorientation in Microgravity, Advances in Space Research . 13(7) . 155-163, 1993. 

(14) Hung, R. J., Lee, C. C., and Leslie, F. W. , Effect of the Baffle on the Asymmetric Gravity-Jitter 
Excited Slosh Waves and Spacecraft Moment and Angular Momentum Fluctuations, Journal of Aerospace 
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(15) Hung, R. J., and Pan, H. L., Asymmetric Slosh Wave Excitation in Liquid-Vapor Interface Under 
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(16) Hung, R. J. , and Lee, C. C. , Effect of the Baffle on On-Orbit Spacecraft Fluid System Angular 
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(17) Hung, R. J., and Pan, H. L. , Differences in Gravity Gradient and Gravity Jitter-Excited Slosh Waves 
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1993. 

(18) Hung, R. J., and Shyu, K. L. , Liquid Hydrogen Shut-Off Geyser Excitation Induced by Sloshing Dynamics 
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II. Basic Characteristics of Gravity Gradient 
and G-Jitter Accelerations 

Any fluid element inside the on-orbit spacecraft fluid system is subject 
to the acceleration that arises from the gravity gradient of the Earth 1-3 . This 
acceleration can be calculated based on the non— inertial frame spacecraft bound 
coordinate. Thus, the coordinate system shall be transformed from ordinary 
inertial frame coordinate to non-inertial coordinate. 

(A) Orbit Motion of Spacecraft 

Let us consider the examples of Gravity Probe— B (GP— B) and Advanced X-Ray 
Astrophysics Facility-Spectroscopy (AXAF-S) spacecrafts, which are the Earth 
satellites orbiting at 650 km altitude directly over the poles, the orbit period, 
r c can be computed from following expression: 


x 0 =2jc 


R 


3/2 

c 




( 2 - 1 ) 


where R K denotes the radius of Earth (*= 6373 km); R*., the radius of the circular 
orbit (- R K + h - 7023 km); h, orbit altitude (« 650 km); and go, Earth gravity 
acceleration (- 9.81 m/s 2 ). For the cases of both GP-B and AXAF-S spacecrafts, 
the orbit period r 0 - 97.6 min, and orbit rate n * - 1.07 x 10“ 3 rad /« • 

As the spacecraft is orbiting around the Earth, the azimuth angle of the 
Earth , toward the location of the spacecraft mass center varies with respect 
to time. At time t — 0, the rolling axis of the spacecraft is aligned with the 
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radial direction of the Earth's center to the spacecraft mass center. Assuming 
the spacecraft rolling axis is linearly turning around 0° to 360° in the orbit 
period, r c , of the spacecraft when the spacecraft is orbiting around the Earth. 
This is particularly true for the case of the GP— B spacecraft. Without the 
spacecraft slew motion, the azimuth angle (^ Eo ) can be defined as 



where r c is the spacecraft orbit period [defined in Equation (2-1)] ; and t is the 
time measured from the instant when the direction of the spacecraft rolling axis 
is aligned with the radial direction of the spacecraft mass center to the center 
of the Earth. 

(B) Slew Motion of Spacecraft 

For the purpose to carry out wide-range observations, some scientific 
spacecraft requires slew motion with respect to the mass center of the 
spacecraft. This is particularly true for the case of the AXAF-S spacecraft. 
For the case of the spacecraft slew motion, azimuth angle, shown in Equation (2- 
2) , shall be modified through the coordinate transformation of slew motion when 
the spacecraft is orbiting around the Earth. 

Let us assume that the slew motion starts with the center located at the 
mass center of the spacecraft. Let us choose cartesian coordinate (x n , y n , z") 
with z"-axis along the axis of the dewar container (see Figure 1). At time t - 
0, the radial vector r c from the spacecraft mass center to the center of the 
Earth lies on the x"-z" plane of the cartesian coordinate chosen (see Figure 1) . 
The azimuth angle defined as the angle between the radial vector r c and the 
z w -axis. Rotation matrices for spinning and/or slew motions along the x”-, y"— 
and z"-axes can be expressed as 
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10 0 


’cosa) y t 0 -sinci> y t 


coso> z t sinco z t 0 

0 coso> x t sinw x t 


0 10 


-sinci> z t coso> z t 0 

0 -sino> x t cosco x t 


sino) y t 0 cos o> y t 


0 0 1 


respectively. Here, w x , w y and denote angular velocity of slew and/or 
spinning motions along the x'-, y'- and z'-axes, respectively. Radial vector r c 
in cartesian coordinate without slew and spinning motion is (see Figure 1) 

- r . . . , ( 2 - 3 ) 

f co = 0 , -cosi|f Ko ] 


For the successive operations of the spacecraft from spinning motion along the 
z"— axis , then slew motion along the y"-axis , and then slew motion along the x w - 
axis, radial vector r c results 


1 

0 

1 


COS G) y fc 

0 

-sinco y t 

0 

coso> x t 

sinu x t 


0 

1 

0 

0 

-sinco x t 

COSQ x t 


sino) y t 

0 

COS(O y t 


’ coso> r t sinco^fc 0 
- -sinco z t coso) r 0 
0 0 1 

In addition to the modification of the azimuth angle made by the spacecraft slew 
motion through the formulation of coordinate transformation, shown in Equations 
(2-3) and (2-4), accelerations are also induced to activate on the fluid mass in 
the dewar container. Accelerations acting on the fluid particle in the dewar 
induced by the slew motion of the spacecraft with the coordinate fixed at the 
spacecraft mass center mass is as follows (see Figure 1): 

R p = Qx(&xRp) + axR p + 2 d XV ( 2 - 5 ) 

► 

where denotes the position vector of the fluid particle in the dewar container 
relative to the body frame of the spacecraft; w, angular velocity of the 


sinilr*, 

0 

-cos^ 


( 2 - 4 ) 
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spacecraft body frame; a, angular acceleration of the spacecraft body frame; and 
v, velocity of the fluid particle relative to the spacecraft body frame. 

As we indicated earlier, let us assume that the slew motion starts with the 
center located at the spacecraft mass center, cartesian coordinate (x' 9 , y ' 9 , 
z'') is chosen with origin located at the spacecraft mass center. Let us also 
assume that x' z ' * plane intersects the center of Earth and the spacecraft mass 
center. In other words, azimuth angle of Earth toward the spacecraft mass center 

lies in the x''-z'' plane. Slew motion is along both the x''- and y''- 

^ 

coordinates. Thus, u a - (cj # x , u> sy , 0) and a m - (a, x , a x y, 0), Rp due to slew motion 
becomes 

„ [ 4"1 [ +a. ry R z +2(J ay v z 

* Ry * ( 2 - 6 ) 

K"J »i.v 1 _J? * ( U « +0> «) + ( a «r R y _a *jr R x) + 2 ( «„V y -U ^,V X ) J ^ 

(C) Coupling for the Accelerations of Spinning and Slew Motion of Spacecraft 
As indicated in this review, with their specified functions of scientific 
observation, the dewar container of the GP-B is spinning with a certain rotating 
rate without slew motion during the normal operation while AXAF-S requires slew 
motion for pointing control to observe point and extended sources of astronomical 
objects without spinning. For some particular reasons required in other 
spacecraft, it might be faced with the situation that both spinning and slew 
motions are needed simultaneously. To encounter this case, the following 
formulations are made to deal with coupling for the accelerations of spinning and 
slew motion of the spacecraft: 

Blew and spinning' 
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^ syftx^ sy^z 

” ^ sx ^ ^ sx^y ~ ^ sy^x ^ "®sx^r”^^sx V x 
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spinning and coupiing 


(2-7) 


where w x and <I) r denote angular velocity and angular acceleration, respectively, 
of spacecraft spinning motion along the z-axis. 

For the case of the GP-B spacecraft, there is no slew motion and the 
spinning is the only acceleration acting on the spacecraft fluid system. 
Acceleration due to spacecraft spinning motion becomes 


R 


■p, spinning- 
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-R o )t-i? cl) -2(0 V 

Iv x ^ y z zy 
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Ry 
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( 2 - 8 ) 


To convert the expression of Equation (2-8) in cartesian coordinate to 
cylindrical coordinate, by using the relationships of (R*, Ry) - (rcostf, rsintf) 
and (v x , v y ) - ( u,. cosfl - u$ sin0, Uj sin0 + u 9 cos$) , Equation (2-8) becomes 


and 


K p, spinning 
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-rcos0o> x -rsin0d) x -2 (i^sinO+i^cosB) o> x 
-rsin0<*>^+rcos0u x +2 (u^osO-i^sinO) <*> x 
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spinning 


(2-9) 
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Accelerations induced by spacecraft spinning motion alone becomes 
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( 2 - 11 ) 


R r ro>|+2u e o r 

= -Rq = =ru r 2u r u r 

R 0 

spinning L rJ spinning spinning 

(D) Gravity Gradient Acceleration 

The gravity gradient acceleration acting on the fluid mass of spacecraft 
can be shown as 

^=-n 2 [ 3 (£ c ^d)£ c -d] ( 2 - 12 ) 

where a gg denotes gravity gradient acceleration vector; the vector (not a unit 
vector) from the fluid element to the spacecraft mass center; r c , a unit vector 
from the spacecraft mass center to the center of the Earth; and n, the orbit rate 
(see Figure 1). 

It is assumed that the gravity gradient exerted on the mass center of the 
spacecraft orbiting around the Earth on its specified orbit is zero. In other 
words, all the gravity acceleration exerted on the spacecraft is nothing but the 
gravity gradient acceleration which is defined in Equation (2—12) . In this 
study, we are interested in investigating how gravity gradient acceleration 
affects the dynamical behavior of cryogenic fluid elements of helium. 

For the convenience of mathematical study, let us describe all the 
parameters Involved in Equation (2-12) in terms of cartesian coordinates. In 
order to match with the computer simulation, mathematical derivation are 
considered in the first quadrant. Figure 1 illustrates the geometrical 
relationship of the parameters shown in Equation (2—12) . 

Let us consider the fluid element of interests, m, located at (r, 0, z) in 
cylindrical coordinates and at (x, y, z) In cartesian coordinates. The origin 
of the two coordinate systems Is located at the center bottom of the dewar tank. 
The slew and/or spinning motions, mentioned earlier, are executed at the 
spacecraft mass center with cartesian coordinate (x" , y" , z"). The spacecraft 
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mass center is located at z - 4* Assume that vector ? c lies in the x-z plane 
of the cartesian coordinate. Vector fi in (x, y, z) coordinate becomes 

cf = [-zcos0, -rsin0, - (z-L c ) ] (2-13) 

Substituting Equations (2-4) and (2-13) in (2-12), non-inertial frame 
expression of gravity gradient acceleration with spinning motion in z-axis 
becomes 


[ &„ m J J3 [ -rain+seCOB (0+u, t) ♦ ( z-L e ) cosf J BiT\+ Bo coB<* B t+rcos& 

J --n 2 -3 [rsin^coB (0+«,t) ♦ (z-L c ) cos «int J& flinw,t>rsin6 

t* ' 3 b ♦(^-X, c )co«t^]co«t Jb ^(z-r c ) 


(2-14 ) 


Similarly, non-inertial frame expressions of gravity gradient acceleration 


with slew motion in y-axis becomes 




3 [-rslntycosd (z-L c ) ] sini|r +rcos0 

^gg ,y 

~-n 2 

rsinG 

a 99. z. 

a law in y-xxia 

-3 [-rsim(rcos0+cosi|r (z-L c ) ] cosiJr + (z-L c ) 


where - V> Eo + «yt. 

Also, non-inertial frame expressions of gravity gradient acceleration with 
slew motion in x-axis becomes 


a 99.x 


BAsimlr^+rcosB 

a 99.y 

= -n 2 
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a gy. z. 

slew in x-axis 

3Acos\Jf B0 coso> x t+ (z-L c ) 


A=-rcos0sini{j Eo +cosi|f Eo [rsinw x tsin0+cosco x t (z-L c ) ] 

The gravity gradient acceleration located at (r, $ , z) can be computed from 
that located at (x, y, z) , from the following relation: 



^39. r 
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(E) G-Jitter Accelerations 

In addition to gravity gradient acceleration acting on the fluid element 
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of on-orbit spacecraft fluid systems, there is another acceleration of g-jitter 
also exerted forces on the fluid systems. 

Among varieties of g-jitter accelerations listed, accelerations induced by 
slew motion of the spacecraft dominate the forces activated on the spacecraft 
fluid systems. In the derivation of acceleration induced by the slew motion of 
the spacecraft, the coordinate system (x" , y" , z") is fixed at the spacecraft 
mass center. A detailed expression of [R*, Ry, R*] # i* w are shown in Equation (2-6) 
of this review. G-jitter acceleration is a summation of acceleration induced by 
slew motion and others, such as atmospheric drag on the spacecraft, spacecraft 
attitude motions arising from machinery vibration, thruster firing, crew motion, 
etc. Thus, g-jitter acceleration can be expressed as 


** = 
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where f is the jitter frequency (Hz) imposed on the fluid systems of the 
spacecraft. The amplitude of the g-jitter fluctuating component, assumed to be 
one-half the mean value, is based upon Russian experience in Mir experiments*. 

III. Non-Inertial Frame Mathematical Formulation of Fundamental Equations 
Dynamical Behavior of fluid elements inside the on— orbit spacecraft fluid 
systems are strongly modified by the gravity gradient and g-jitter accelerations. 
In order to accommodate the effect of gravity gradient acceleration on the on- 
orbit fluid motion, one has to consider non-inertial frame spacecraft bound 
coordinate rather that adopting inertial frame coordinate used in ordinary fluid 
mechanics formulation. 

Experiment made by Andronikashvili 23,2 * for rotating helium II shows that 
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it is necessary to exceed a critical velocity for the interaction between the 
normal and superfluid components to establish entire bucket in rotation 23 . 
Based on this illustration shown earlier, the problem under consideration with 
normal and superfluid helium, which can produce great enough vortex line tangles 
to assure an interaction between the normal and superfluid components, suggests 
an adoption of Newtonian fluids in terms of corresponding temperature with proper 
concentration ratio of normal and superfluid helium in this review 30 * 31 . 

Consider a closed circular cylindrical dewar of radius, a, with height, L, 
which is partially filled with cryogenic liquid helium, and the ullage is filled 
with a helium vapor. Angular velocity of rotating cylinder is u>. Density and 
viscosity of liquid helium and helium vapor are p L , p v , and respectively. 
Let us use cylindrical coordinates (r, 6, z) , with corresponding velocity 

components (u, v, w) , and corresponding residual gravity acceleration, such as 
gravity gradient components (a ggr , a gg 0t a gg z ) and g-jitter components (a gj r , 
a ff j 0 , a g j t ) . In the derivation of the governing equations, accelerations induced 
by the spinning motion of the spacecraft is included in the formulation. The 
rest of the acceleration such as slew motion, atmospheric drag on the spacecraft, 
spacecraft attitude motions arising from machinery vibrations, thruster firing, 
and others, are included in the jitter acceleration, shown in Equation (2-18). 
The governing equations for non— inertial frame of spacecraft bound coordinates 
can be shown as follows : 

(A) Continuity Equation 


I JL 

r dr 


(ru) + 


1 dv ^dw 
r dQ dz 


=0 


(3-1) 


(B) Momentum Equations 
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where 


p 'Ir-f + Tw-i> =-H * p <’-<> 
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r dr dr r 2 36 2 3z 2 
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In these formulations, 2cj 2 v and 2io z u are the Coriolis acceleration, rw z 2 is 
the centrifugal acceleration, and rcj r is the angular acceleration induced by the 
spinning motion of the spacecraft. 

IV. Initial and Boundary Conditions of 
Spacecraft Fluid System in Microgravity Environment 
Governing equations of the fluid motion in on— orbit spacecraft fluid 
systems for non-inertial frame spacecraft bound coordinates have been illustrated 
in Equations (3-1) to (3-5). These equations shall be combined with the 

characteristics of gravity gradient and g-jitter accelerations as that formulated 
in Equations (2-1) to (2-18) . Initial and boundary conditions shall be 
introduced to accommodate solving fluid motion in on-orbit spacecraft fluid 
system for non-inertial frame coordinate 7 ' 9,32 ” 34 . 

There are several ways to compute the dynamical behavior of liquid— vapor 
interface. One way is explicitly tracking the dynamics of interface fluctuations 
through the computation of the radius of curvature variations along the 
interface 35 . The other method is to implicitly capture the interface using a 
level set approach, and the interface will be identified as the zero level set 
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of a smooth function 36 . In this review, the method of explicitly tracking the 
interface fluctuations is described 35 . 

Let the profile of the interface between gaseous and liquid fluids be given 

by: 

x] ( t, r, 0, z) = 0 (4-1) 

The initial condition of the profile of the interface between gaseous and 
liquid fluids at t - t c is assigned explicitly, and is given by: 

n (C = t Q , r, 0, z) - 0 (4-2) 

A set of boundary conditions has to be supplied for solving the equations. 
These initial interface profiles used in this study have been given explicitly 
through the steady state computations made by Hung and Leslie 20 and Hung et al 7 “ 9 
which were checked by the experiments carried out by Leslie 21 . These boundary 
conditions are as follows: 

(1) Along the container wall, the following three boundary conditions apply: 

(a) Interface between solid and liquid fluid: No-penetration and no— slip 
conditions assure that both the tangential and the normal components 
of the liquid velocity along the solid walls will vanish. 

(b) Interface between solid and gaseous (vapor) fluid: Similar no- 

penetration and no-slip conditions as that shown for interface 
between solid and liquid fluid will apply. 

(c) At the location of solid-liquid-gaseous (vapor) three phases 

interface: No-penetration, but not no-slip condition apply. This 

will assure that normal components of liquid and vapor velocities 
along the solid walls will vanish, and allow a slipping flow of 
liquid and vapor fluids along the solid wall at three phase interface 
location. . The velocity of slipping flow at this location is governed 
by the adhesive forces between fluids (liquid and gaseous) and solid 
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walls. Also, at this location of three phase interface, a constant 
contact angle is present in which the behaviors of wet or dry 
contacts are determined by Coulomb interaction between the fluids 
(liquid and vapor) and the surface phenomena (material and roughness) 
of solid walls 35 ' 39 . 

(2) Along the interface between the liquid and gaseous fluids, the 
following two conditions apply: 

(a) Kinematic surface boundary condition: The liquid (or gaseous) 

surface moves with the liquid (or gas) which implies 

^n=0,or 

Dt 

41 +u-^- + — 4S - =0 on *1 ( t=t «» r,0, z) (4-3) 

dt dr r dd dz 

(b) Interface stress condition: Across the liquid-vapor interface, the 

stress must be continuous. Based on Landau and Lifshitz 40 , the 

stress across the liquid-vapor interface can be expressed as 

where R x and R 2 are the radii of curvatures of two major axes at the 
point of interests on the surface of the liquid vapor interface. 

The expressions of radii of curvatures R x and R 2 in cylindrical coordinates 
from differential geometry can be shown as 40 ' 42 
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and the configuration of the liquid-vapor interface is z - H(t, r f 0). Here, in 
Equations (4-4) and (4-5) 




= M 


du 1 

dx., 


duj 

dx ± 


du * 6 


d u u 

OXy 


ij 


is the viscous stress tensor; the viscous coefficient of the first kind; £, 
the viscous coefficient of the second kind; P, the pressure; a, the surface 
tension of the liquid-vapor interface; and n^ , the unit vector normal to the 
interface; and 6 ^ , the Kronecker delta function. Also, subscripts G and L 
denote the conditions at gaseous and liquid fluids, respectively, across the 
liquid-vapor interface . 

The fluid stresses across the liquid— vapor interface can be decomposed to 
the components normal (n if a unit vector) and tangential (t A , a unit vector) to 
the interface. For the component tangential to the interface, one can take a dot 
product of a unit vector tangential to the interface, t lf to Equations (4—4) and 
(4—5) , which leads to 


[(^ ij t i n j )] L =[x lj t i n :j )] G (4-6) 

since u i t i - 0. 

For the component normal to the interface, one can also take a dot product 
of a unit vector normal to the interface, r^, to equations (4-4) and (4-5), which 
leads to 


p a - Pl - " 


( X lj n i n j) J 



0 r^ r ( 
'Sr~D~ ) 



‘-s’ 1 


(4-7) 


For components normal to the interface along the (r, 0, z) directions in 
cylindrical coordinates can be obtained by taking dot products of n*, n B , n* 
separately to Equations (4-4) and (4—5), which are expressed as 
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where (i^ , n fl , n*) is the unit vector normal to the interface in cylindrical 
coordinates (r, 0, z) . 

For the purpose of solving bubble dynamic problems in microgravity, one 
must solve the governing equations shown in Equations (3—1) to (3-5) , accompanied 
by a set of initial and boundary conditions, shown in Equations (4-1) to (4-8). 
The computational algorithm applicable to cryogenic fluid management under 
microgravity is also given earlier 32 " 37 . Summarized computational algorithm are 
illustrated in Figure 2. Most of the illustrations shown in this review adopt 
a full scale GP-B spacecraft dewar tank with an inner radius of 12 cm and an 
outer radius of 68 cm and a height of 145 cm. The dewar tank is partially filled 
with cryogenic helium and the rest of the ullage is filled with helium vapor. 
The temperature of cryogenic helium is 1.8 K. In this review the following data 
were used: liquid helium density - 0.146 g/cm 3 , helium vapor density - 0.00147 
g/cm 3 , fluid pressure - 1.66 x 10* dyne/cm 2 , surface tension coefficient at the 
interface between liquid helium and helium vapor - 0.353 dyne/cm, liquid helium 
viscosity coefficient - 9.61 x 10" 5 cm 2 /s, and contact angle - 0°. 

A staggered grid for the velocity components is used in this computer 
program. MAC (marker-and-cell) method* 1 of studying fluid flows along a free 
surface is adopted. VOF (volume of fluid) method is used to track interface by 
solving finite difference equations numerically. Figure 2 describes computation- 
al algorithm and numerical method to track liquid— vapor Interface profiles 
numerically. The approximate flow velocity is calculated from the explicit 
approximation of momentum equations based on the results from the previous time 
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step. Computation of pressure and velocity at the new time step are, thus, 
obtained from iteratively solving the pressure equation through conjugate 
residual technique 41 " 44 . The configuration of liquid-vapor interface adjusted by 
the surface tension effect at the new time step are then finally obtained. The 
time step during this computation is automatically adjusted through the 
fulfillment of the stability criteria of computed grid size. Convergence 
criterion of the iteration of pressure equation is based on the computed velocity 
at each cell which satisfy continuity equation with the errors no more than 10 -5 
of the velocity difference 45,46 . As for the volume conservation of liquid, a 
deviation of less than 1 % error of volume is always guaranteed before a move to 
the next time step. 

V. Mathematical Formulations of Slosh Reaction Forces and Torques Due to 

Bubble Deformations 

Sloshing dynamics of fluid system will introduce bubble deformations and 
oscillations which can be characterized by the types of externally applied forces 
in the microgravity environment. These disturbed bubble through the deformations 
and oscillations also produce slosh reaction feedback forces and torques acting 
back on the fluid filled container. Computation of slosh reaction forces and 
torques due to bubble deformation is essential for the development of spacecraft 
orbital, attitude and guidance control techniques. Mathematical formulation is 
derived to proceed the computation of slosh reaction forces and torque 
fluctuations due to bubble deformation 47-50 . 

For the purpose of considering large amplitude slosh wave activated slosh 
reaction forces exerted on the solid walls of the dewar, the fluid stresses are 
decomposed into the tangential and normal components acting on the walls which 
can be expressed as follows: 
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where n t denotes the tangential component of fluid stresses; 1^, the normal 
component of fluid stresses; P, the thermodynamic pressure; fluid velocity 
in a direction; i ai unit vector tangential to the wall; n^, unit vector normal 
to the wall; p, the molecular viscosity coefficient of fluid; and 6 a ^ t the 
Kronecker delta function. Subscripts a and p imply the directions of flow 
fields . 

The stress distribution shown in Equations (5—1) and (5-2) can be 
integrated with respect to area and obtain the tangential and normal forces 
acting on the container wall. Detailed mathematical derivations of slosh 
reaction forces acting on the walls have been fully illustrated 32,35,46,53-56 , and 
will not be repeated in this review. 

In order to accommodate the spacecraft dynamics of pitching, yawing and 
rolling, cylindrical coordinates of the rotating container is transformed into 
cartesian coordinates based on (x,y,z) - (rcos0, rsin0,z) with corresponding 
velocity components (v x , v y , v z ) - (ucosfl - vsintf, usintf + vcostf , w) . For the 
case of the GP-B spacecraft, the axis of rotating is always fixed at the point 
of proof mass which is located at the geometrical center of the dewar at (x c , y c , 
z c ) - (0, 0, Lc) , where - 1/2 L and L is the height of the dewar. For the 
case of the AXAF-S spacecraft which is located at (x Cf y c , z c ) - (0, 0, 1^) where 
Lc is the height of the axis of slew motion (see Figure 1) . 

In addition to the slosh reaction forces (F x , F y , F x ) along (x, y, z) 
coordinates, slosh reaction torques (M*, My, M*) acting on the dewar container 
also have been derived mathematically 32,35,46,53-56 , and will not be repeated in 
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this review. Relationships of slosh reaction forces (F x , F y , F r ) , slosh reaction 
torques My, M x ) and slosh reaction moment arms (L*, Ly, L^) are as follows: 


\M X ] 

fV^-(L,-L c )F y j 


M y 

= (L M -L e )F x -L x F t 

(5-3) 
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By using the relation of F a L* - 0, or F x L* + F y Ly + F x (L* - Lc) - 0, moment arms 
of slosh reaction torques induced by the slosh wave excitation can be computed 


from the following relations: 


A, 1 - \FJt,-FM 

L y - * -i F A-FA 
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(5-4) 


The problems illustrated can be applied to bubble deformations activated 
by various forms of accelerations. Following examples are given to illustrate 
the bubble deformations driven by (a) centrifugal forces — bubble spin-up and 
spin-down motions, (b) lateral and axial impulses, (c) slew motion accompanied 


accelerations, and (d) combined gravity gradient and g-jitter accelerations. In 
these applications, first of all the time sequence animation of bubble 
deformations, driven by the various forms of accelerations mentioned earlier, 
will be illustrated based on the interests of the fluid dynamics point of view. 


Bubble mass center fluctuations together with slosh reaction forces and torques 
Induced due to bubble deformations will also be described based on the timely 
interests from spacecraft orbital, attitude, and guidance control technology 
development points of view. 


VI. Bubble Spin— Up From Rest With and Without a Completion of Wrapping 

Bubble is Initially sitting at the location between outer and inner dewar 
wells without rotation. During the process of bubble deformation from spherical- 
shape to the final shape of bubble wrapping around the inner well of the rotating 
dewar, the critical point Is that the rotating speed of the dewar shall be great 
enough to create a sufficiently large centrifugal force to overcome surface 
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tension of bubble such that the bubble can stretch, elongate, and finally wrap 
around the dewar well. Dimensionless parameters characterized by the Weber 
number (We) i s 37 .50.5i 

We = P Q2r3 = centrifugal force 

o surface tension force 

Here, p, r, and a denote fluid density, angular velocity, radius of the 
container, and surface tension of bubble at the interface between liquid and 
vapor helium. In other words, We number also characterizes the bubble wrapping 
around the dewar inner well during the spin— up motion. A greater value of We 
number is in favor of a complete wrapping of the bubble around the dewar inner 
well 32 . 

To show an example of bubble spin-up with a complete wrapping for 95% 
liquid-fill level, the dewar container is spinning up from w - 0 to 0.4 rpm 
linearly in 2 s and keeps constant angular velocity thereafter. Figure 3(A) 
shows the time sequence evolution of the three dimensional dynamical behavior of 
bubble evolution through the processes of elongation, wrapping around the Inner 
well of the dewar, and finally a completion of the bubble wrapping around the 
inner well during the spin-up rotating of dewar. It shows that the bubble 
gradually deforms from spherical shape to doughnut shape which wraps around the 
inner column of the dewar, and becomes quasi-symmetric profile at time t - 8700 
s and thereafter. 

For the purpose to check the results of numerical simulation, dimensionless 
time scale of E k “ 1/2 is used to estimate time required for the development of 
viscous boundary layer 33 . Parameter E* is Ekman number which is defined by 33 


= 


0) (A R) 


( 6 - 1 ) 


Here v stands for kinematic viscosity, and AR - (r„ - r A )/2 where r e and r A denote 
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outer and Inner radii, respectively, of the dewar. In reality, Ekman number is 
nothing but the reverse of Reynolds number for rotating system (R*” 1 ) . Estimated 
time required for the fluid to attain the angular velocity of new steady state 
is on the order of E k ~ 1/2 /u during spin-up motion 53 . To show examples, estimated 
time required for cryogenic helium to reach steady state rotating speed of 0.4 
rpm spin-up from rest is on the order of 10 4 s while that of the rotating speed 
of 0.2 rpm spin-up from rest is on the order of 1.4 X 10 4 s. This order of 
magnitude estimation of time required to attain steady state is agreeable with 
the numerical simulation shown in this review. 

Weber number also characterizes the time required for bubble wrapping to 
reach equilibrium shapes. In other words, the higher the Weber number, the 
shorter the time required for the bubble to reach equilibrium shapes . It shows 
that the equilibrium bubble profiles with liquid fill level of 95% for bubble 
wrapping can be resulted with (a) We-56.65, w-0.2 rpm, and time required to reach 
equilibrium at t-13,100 s; (b) We-14.18, c^-0 . 1 rpm, and time at t-34,200 s; and 
(c) We-3.54, a>-0.05 rpm, and time at t-87,000 s. Degrees of bubble wrapping with 
different Weber numbers for liquid fill level of 95% are illustrated in Figure 
3(B). It shows that higher Weber number (higher rotating speed) suggests a 
higher degree of bubble wrapping. This figure shows that the minimum rotating 
speed for a completion of a full wrapping is 0.2 rpm with liquid fill level of 
95%. In other words, bubble wrapping with a rotating speed less than 0.2 rpm 
will be ending at an incomplete bubble wrapping. Figure 3(B) also shows various 
angles of wrapping with liquid fill level of 95% for an incomplete bubble 
wrapping corresponding to different rotating speeds with a><0.2 rpm. Angles of 
incomplete wrapping are 53°, 124°, and 178° for bubble wrapping with rotating 
speeds corresponding to 0.15, 0.1, and 0.05 rpm, respectively, with liquid fill 
level of 95%. 

With reference to Figure 3(A) for bubble spin-up from rest to rotating 
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speed of 0.4 rpm with liquid fill level of 95%, Figure 4(A) shows the time 
evolution of bubble mass center fluctuations during the process of bubble 
wrapping around the dewar well. The figure shows the bubble mass center time 
evolution x c , y c and z c along the x, y and z axes, respectively, in non-inertial 
bubble frame. This figure provides the following conclusions for bubble mass 
center fluctuations: (a) At time t - 0, the bubble with the spherical-shape was 
sitting at x c - -38.2 cm, while y c and z c are basically located at a central 
position. (b) As the dewar container starts to spin-up, the bubble mass center 
along the x— axis gradually moves toward a central location while the bubble mass 
center along the y-axis shifts slightly from a central location with small 
amplitude oscillations and bubble mass center along the z— axis with negligible 
amplitude of oscillation. (c) The major driving forces for the oscillations of 
bubble mass fluctuations along the x- and y-axes are due to the centrifugal and 
Coriolis forces while none of these forces apply to the z-axis. (d) At the very 
beginning, x c oscillates from - 38.2 cm toward the central position while y c 
fluctuates from the central position toward negative y position, (e) Oscillation 
periods and amplitudes for bubble mass center fluctuations along the y-axis are 
such that the initial mode with a shorter periods (2000 s) and a greater 
amplitude (5 cm) and subsequent second mode with longer period (5000 s) and a 
smaller amplitude (3 cm). In other words, the bubble mass center oscillations 
along the y-axis change from a shorter period to a longer period, and from a 
greater amplitude to a smaller amplitude, and then gradually dies out. (f) 
Similar characteristics of the variation of a period and an amplitude for bubble 
mass center fluctuations and oscillations along the x-axis are shown which is 
similar to that of the figures along the y— axis . 

Furthermore, one can study the time evolution of bubble mass center 
fluctuations with and without a completion of wrapping around the dewar well 
driven by the spin-up motion. Figures 4(B), 4(C), and 4(D) show bubble mass 


28 


center evolution with rotating speeds of 0.2, 0.15, and 0.05 rpm, respectively. 
These figures provide the following conclusions for bubble mass center 
fluctuations in addition to the conclusions given for Figure 3(A): (a) The 
bubble mass center oscillations along the x-axis finally move to a central 
location for the container rotating speeds w>0.2 rpm while bubble mass center 
oscillations are ending at some equilibrium positions near but not at the central 
position for the container rotating speeds cj< 0.2 rpm. (b) A greater gap between 
bubble equilibrium locations and central location along the x-axis is shown for 
the lower container rotating speeds or lower We number in comparison with cj— 0.2 
rpm and We— 56.65 for minimum values required for the bubble to accomplish a 
complete wrapping process for the container liquid fill level of 95%. 

Bubble spin-down from a complete wrapping with steady state is also 
investigated. At time t— 0, the dewar container is spinning with 0.2 rpm with a 
full bubble wrapping around the dewar well for liquid-fill level of 95%. Figure 
5(A) shows the time sequence evolution of the three dimensional bubble 
configuration and select times for spin— down motion of dewar with rotating speeds 
from 0.2 to 0.15 rpm. Bubble behaves dynamical evolution through the processes 
of deformations in terms of shrinking in the direction longitudinal to, and 
stretching in the direction transverse to the rotating axis, jointly, without 
unwrapping around the dewar well during spin-down motion. 

Further decreasing of rotating speed during the spin-down motion than that 
illustrated in Figure 5(A) eventually causes the motion of bubble unwrapping 
around the dewar well for dewar with liquid fill level of 95%. Figure 5(B) shows 
the time sequence evolution of the three dimensional configuration at some select 
times for spin-down motion of dewar with rotating speeds from 0.2 to 0.05 rpm 
when bubble starts to unwrap around dewar well with liquid fill level of 95%. 
Figure 5(B) also clearly shows that the symmetric rotating bubble wrapping around 
the dewar well starts to create an asymmetric ripple along the surface of the 
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bubble at t-4,650 s. These ripple perturbations also quickly enhance and 
propagate. Finally the weakening centripetal forces resulted during the spin- 
down motion are unable to maintain the bubble wrapping around the dewar well. 
Thus, asymmetry in the imbalance liquid-vapor distribution is then starting to 
create . 

There are some basic differences between the realistic practice and the 
theoretical computation for the problems of bubble spin— down motion in rotating 
dewar. In the theoretical computation, a symmetrical profile of distribution 
with respect to the rotating axis for centrifugal and Coriolis forces together 
with angular acceleration, which are all nothing but the function of angular 
velocity of rotating speed, shall be assumed during the time period of spin— down 
motion. In the realistic practice, a pronounced asymmetric force will be induced 
in these weakening forces of spin-down motion and initiates a cracking at this 
particular location of a symmetric force activating on the wrapped bubble which 
is the major cause of crack propapation leading to the final unwrapping of the 
rotating bubble around the dewar well. Thus, the symmetrical wrapping bubble 
will end up in unwrapping status during the spin-down motion. In the theoretical 
computation, one can never know where and when an asymmetric force shall be 
initiated unless this force is added into the computation artificially and 
arbitrarily. This Is because there is no mechanism In the theoretical modelling 
to initiate such an asymmetric force to unwrap the rotating bubble to balance the 
weakening centrifugal forces during the spin-down motion. In other words, the 
basic differences between the realistic practice and computation are (a) the 
bubble will be unwrapped and still maintain In one piece or even two pieces 
encircling around the rotating dewar well depending upon the strength of surface 
tension which has a strong intention to maintain a single piece of rotating 
bubble for the case of realistic practice, and (b) the bubble will be broken into 
pieces because of the axial symmetric forces acting on the bubble with weakening 
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magnitude of centrifugal force which is unable to ensure bubble wrapping around 
the rotating dewar well for the case of theoretical computation 51 . 

The Weber numbers considered in the cases of bubble evolution, shown in 
Figures 5(A) and 5(B) are 31.87 and 3.54, respectively. It shows that the Weber 
number characterizes the possibility of the bubble with and without unwrapping 
around the dewar well. Results indicate that the rotating bubble can maintain 
its symmetric configurations without unwrapping around the dewar well if the 
Weber number is greater than 3.54 with dewar liquid fill level of 95%. 
Otherwise, the centripetal forces are too weak to maintain the rotating bubble 
without unwrapping around the dewar well. 

A comparison of rotating dewar with liquid fill levels of 85 and 95% shows 
the results with the following conclusions: (a) The higher the liquid fill level 
of the rotating dewar, the higher the rotating speed required for a bubble to 
complete a full wrapping around the dewar well during spin— up motion (0.2 and 
0.08 rpm rotating speeds are required for bubble to complete a full wrapping 
around dewar well with liquid fill levels of 95 and 85%, respectively.) (b) The 
higher the liquid fill level of the rotating dewar, the higher the rotating speed 
required for the bubble to maintain unwrapping around the dewar well during the 
spin— down motion (0.05 and 0.02 rpm rotating speeds are required to maintain a 
full wrapping bubble without unwrapping around the dewar well with liquid fill 
level of 95 and 85%, respectively). (c) The higher the liquid fill level of 
rotating dewar, the shorter the time required for the rotating bubble to 
accomplish a full wrapping around the dewar well (t-13,100 and 30,400 s are 
required for the rotating bubble to complete a full wrapping around the dewar 
well with liquid fill levels of 95 and 85%, respectively). (d) The higher the 
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liquid fill level of the rotating dewar, the shorter the time required for the 
rotating bubble during spin-down motion with full wrapping to become unwrapped 
from the minimum rotating speed required to accomplish a full wrapping during 
spin— up motion (t-4,650 and 10,600 s are required for the rotating bubble with 
full wrapping to become unwrapped from the minimum rotating speed required to 
accomplish a full wrapping around the dewar well with liquid fill levels of 95 
and 85%, respectively, during the spin-down motion). 

Comparison between minimum rotating speeds required for bubble to complete 
a full wrapping during spin-up motion, and minimum rotating speeds required for 
bubble to maintain a full wrapping without unwrapping around the dewar well 
during spin-down motion is quite different. It shows that a minimum speed to 
accomplish a full wrapping during spin-up motion, and a minimum speed to maintain 
a full wrapping without unwrapping during spin— down motion are 0.2 and 0.05 rpm, 
respectively, with liquid fill level of 95%, while the minimum rotating speeds 
for similar spin-up and spin-down motions are 0.08 and 0.02 rpm, respectively, 
with liquid fill level of 85%. These discrepancy in two figures between spin-up 
and spin-down motions are due to the fact that it needs additional energy to 
change any status of the most stable shape with minimum surface area required by 
the condition of minimum Gibbs free energy. The situation of the existence of 
energy barrier between these two motions which are very much similar to the 
motion of electron such as the thermionic— diode generator of the space— charge- 
barrier energy between cathode and anode. 

VII. Rotating Bubble Subject to Lateral and Axial Impulses 

In this review, the response of rotating bubble to impulsive disturbances 
is studied. In order to reduce the degree of asymmetry in the liquid-vapor 
distribution and damp its associated disturbances, a number of baffle are 
installed inside the dewar 50 • 54 * 55 . The difference in response on dewar with and 
without the baffle is also investigated. 
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Figures 6(A) f 6(B) , and 6(C) show three-dimensional configuration of baffle 
installation, the distribution of grid points for the dewar tank with baffle 
boards in the radial-axial plane and radial-circumferential plane, respectively, 
in cylindrical coordinates 56 " 60 . The baffle are installed with three parallel 
plates perpendicular to the rotating axis and four radial plates aligned with the 
rotating axis. The inner radius, outer radius and thickness of the three 
parallel baffle chosen in this review are 60, 77.48 and 0.1016 cm, respectively. 
These three baffles are installed at locations z lt z 2 , z 3 of 38.74, 80.94, and 
84.4 cm, respectively. Four radial baffle-boards are located at 90° apart with 
the same inner and outer diameters and thickness as that of the three parallel 
baffles . 

An example is given to illustrate sloshing dynamics in response to lateral 
impulse during the time period of guidance and/or attitude controls 63 " 66 . Dewar 
is spinning with a rotating rate of 0.1 rpm during its normal operation. In this 
review, with the following form of lateral impulse is assumed: 

a - [a,,, ay, a z ] - [1, 0, 0] 10"^ in cartesian coordinates 
a - [a r , a 0 , a z ] - [cos0, -sinfl, 0] 10"^ in cylindrical coordinates for 
0<t<10" 2 s 

and a - [0. 0. 0] for t>10" 2 s 

where go (9.81m/s 2 )is the Earth gravitational acceleration. 

The time evolution of the bubble fluctuations without a baffle is first 
investigated. In order to illustrate that the bubble were experiencing finite 
deformation, Figure 27(A) shows the time evolution between time from 225 to 427 
s for three-dimensional bubble oscillations without baffle. These three 
subfigures show a transition from left-side bubble swelling, equilibrium and 
right-side bubble swelling during these periods of time. 

In this review, the time evolution of bubble oscillations for the container 
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with baffle due to the sloshing dynamics driven by the same impulse as that 
actuating on the dewar without baffle has been investigated. Figure 27(B) also 
illustrates the time evolution between time from 225 to 427 s for three- 
dimensional bubble oscillations for rotating dewar with baffle to show that 
bubble was experiencing finite deformation with a transition from left-side 
bubble swelling, equilibrium, to right-side bubble swelling. 54,55 

Another example is given to illustrate sloshing dynamics in response to 
axial impulse during the time period of guidance and/or attitude controls. 60,61 
Dewar container is also spinning with a rotating rate of 0.1 rpm during its 
normal operation. In this review, the container bound coordinates with the 
following form of axial impulse are assumed: 

a - [a*, ay, a 2 ] - [0, 0, 1] 10" 2 g o in cartesian coordinates 
a - [a r , a 0t a 2 ] - [0, 0, 1] lO^gQ in cylindrical coordinates for 0<t<10~ 2 s 
and a - [0. 0. 0] for t>10“ 2 s 

Because of page limitation, time sequence evolution of bubble deformations 
activated by the axial impulse without and with baffle will not be illustrated 
in this review. Results of computer simulation show that the bubble is inflated 
at the lower half of the bubble as liquid is pushed upward along the axial 
direction at the very beginning; and then the bubble Is inflated at the upper 
half of the bubble as the liquid is pushed downward along the axial direction 
during the up and down oscillations activated by axial impulse. 

Bubble mass center fluctuations driven by lateral and axial impulses have 
also been studied. Figures 8(A) and 8(B) show the time evolution of bubble mass 
center at (x c , y c , z c ) driven by the lateral and axial impulses, respectively. 
Figure 8(A) shows the time evolution of the growth and decay of the bubble mass 
center fluctuations with solid lines for the rotating dewar with baffle in 
response to the impulse. Figures 8(A-a) , 8(A-b), and 8(A-c) show the time 
evolution of the locations of the fluid mass center fluctuations In x c , y c , and 
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fluctuating locations in radial axis of cylindrical coordinates [ - ( x c 2 + 
y c 2 ) 1/2 ], respectively, for the container with baffle. Similar to the container 
without baffle, the effect of Coriolis forces which carry over the fluctuations 
from x- to y-axes also applies to this case. This means that the fluid mass 
center fluctuations are primarily in the x- and y— axes and not in the z— axis. 
The locations of maximum absolute values of the fluid mass center fluctuations 
in the radial direction are Max(l xj , I yj )-(2.24, 1.67) cm; while the 
fluctuating values of the maximum fluid mass center disturbances in the radial 
coordinates are Max(Ax c , Ay c )~(2.64, 2.52) cm. Similarly, it also shows Ax c >Ay c 
for fluctuating fluid mass center disturbances in rotating dewar with baffle in 
response to the impulse along the x-axis. 

A comparison between the container with and without baffle also illustrates 
some peculiar behavior of cryogenic helium fluids with temperature below A— point 
(2.17 K) in which helium demonstrates a number of remarkable properties of 
superfluidity such as extremely low viscous and surface tension coefficients in 
response to impulse in microgravity environment. It can be concluded as follows: 
(a) Impulse started to actuate on the fluid system at time t - 0 with peak value 
and end at t - l(T 2 s for both the container with and without baffle. However, 
the fluid mass center fluctuations reacted gradually in response to the impulse 
with a maximum amplitude of fluid mass center fluctuations in the x-axis about 
the same time as that between the container with and without baffle at time t - 
80 s, there is a 15 s delay of maximum amplitude fluid mass center fluctuations 
in the y-axis for the container with baffle than that for the container without 
baffle. (b) Displacement of fluid mass center was sinusoidally shifting 
continuously for both the container with and without baffle even long after the 
ending of the impulsive force. This displacement of fluid mass center in the 
radial direction for both cases reached the maximum value at 82 s after the 
ending of impulsive force. Then the fluid mass center started to oscillate back 
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and forth with a longer period of oscillations for the container with baffle than 
that for container without baffle due to the effects of surface tension forces 
and viscous forces for the container with baffle, in addition to the reflection 
and bouncing of flows from the wall. (c) The oscillations of fluid mass center 
fluctuations continued for a long period of time even after the ending of 
impulsive forces due to extremely low viscosity of helium II fluid which 
eventually exponentially decays the mass center fluctuations. (d) As for the 
effect of viscous force between the liquid and solid interface in the container 
with baffle, it contributes a greater damping effect of mass center fluctuations. 
A comparison of the container with and without baffle shows that the container 
with baffle can reduce 13% of the maximum value of fluid mass center dislocation 
and decrease 25% of the maximum amplitude of fluid mass center fluctuations than 
that for the container without baffle. 

Figure 8(B) shows the time evolution of the growth and decay of the fluid 
mass center fluctuations with broken lines in the rotating dewar without baffle 
in response to the axial impulse. The fluid mass center at time t - 0 is located 
at (x c , y c , z c ) - (0, 0, 80.9) cm. The highest amplitude fluctuations of fluid 
mass center is located at Max (x c , y c , z c ) - (0, 0, 86.3) cm at time t - 138 s, 
while the lowest amplitude fluctuations of fluid mass center is located at Min 
(x c , y c , z c ) - (0, 0, 76.6) cm at time t - 406 s for dewar without baffle. In 
this simulation, it shows that fluid mass center fluctuations reach to its 
maximum values at 138 s after the ending of the impulse actuating on the rotating 
dewar. The wave period of fluid mass center fluctuation is about 500 s for dewar 
without baffle . 

Figure 8(B) also shows the time evolution of the growth and decay of the 
fluid mass center fluctuations with solid lines in the rotating dewar with baffle 
in response to the impulse. The highest amplitude fluctuations of fluid mass 
center is located at Max (x c , y c , z c ) - (0, 0, 84.2) cm at time t - 161 s, while 
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the lowest amplitude fluctuations of fluid mass center is located at Min (x c , y c , 
z c ) — (0, 0, 80) cm at time t - 863 s for dewar with baffle. In this simulation, 
it shows that fluid mass center fluctuations reach to its maximum at t - 161 s 
after the ending of the impulsive thrust actuating on the rotating dewar. The 
wave period of fluid mass center fluctuation is about 1000 s. 

A comparison between the container with and without baffle In response to 
impulse in microgravity environment it can be concluded as follows: (a) Impulse 

started to actuate the fluid system at time t - 0 with peak value and end at t 
- 10~ 2 s for both the container with and without baffle. However, the fluid mass 
center fluctuations reacts gradually in response to the impulse with a maximum 
amplitude of fluid mass center fluctuations in the z-axis at t - 138 s for dewar 
without baffle, and at t - 161 s for dewar with baffle, a 23 s delay for dewar 
with baffle to that for dewar without baffle, (b) The oscillations of fluid mass 
center fluctuations continue for a long period of time even after the ending of 
impulse due to extremely low viscosity of helium II fluid which eventually decays 
the mass center fluctuations. (c) As for the effect of drag force between the 
liquid and solid interface in the container with baffle, the installation 
contributes greatly in the increasing of wave period oscillation for the dewar 
with baffle, which shifts the wave period from about 500 s for dewar without 
baffle to about 1000 s for the dewar with baffle. (d) The baffle also 
contributes to the reduction of fluid mass center fluctuations. Results show 
that the container with baffle can reduce 40% of the maximum amplitude of fluid 
mass center fluctuations and decrease 70% of the minimum amplitude of fluid mass 
center fluctuations than that for the container without baffle. 

VIII. Bubble Deformations Subject to Gravity or G— Jitter Acceleration 

Associated with Slew Motion 

By using the mathematical formulations illustrated in Sections II and III, 
one can numerically simulate bubble sloshing dynamics associated with slew 
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motions depending upon the specific scientific missions assigned to the 
spacecraft. An example is given to illustrate sloshing dynamics associated with 
slew motion for the AXAF-S spacecraft. 

(VIII-A) Bubble Deformation Driven by Gravity Gradient Acceleration 

Associated with Slew Motion 

As indicated earlier, the AXAF-S spacecraft is capable of observing point 
and extended sources of active galactic nuclei, clusters of galaxies, supernova 
remnants, x-ray binaries, etc., through the spacecraft slew motion of pointing 
control. Assuming that slew motion is along the y"— axis (see Figure 1), gravity 
gradient acceleration associated with slew motion can be computed from the non- 
inertial frame expressions of Equations (2-12) and (2-17). It is assumed that 
the slew motion operates at 90° in 600 s. 

In this example, spacecraft sloshing dynamics driven by gravity gradient 
acceleration associated with slew motion in the y"-axis, shown in Figure 1, have 
been investigated. As the orbital period of AXAF-S spacecraft is 97.6 min and 
period of 90° slew motion in the y"— axis is 600 s, the component of gravity 
gradient acceleration along the (x, y, z) directions acted on the fluid mass 
located at (r, 0, z) - (12 cm, w/2, 3 cm) is shown in Figure 9(A). This figure 
shows that the magnitude of gravity gradient acceleration is on the order of 10~ 7 
go for AXAL-S dewar on its operation orbit. The distance from the spacecraft 
mass center to the bottom of the dewar, L^, shown in Figure 1, is 257.8 cm. 

The equilibrium shape of the liquid-vapor interface for a dewar with 70% 
liquid-fill level under a residual gravity environment below 10“ 7 go, Is a 
sphere. Figure 10(A) shows the time sequence evolution of the three-dimensional 
dynamical behavior of the bubble oscillations driven by gravity gradient 
acceleration associated with slew motion. The figures of liquid-vapor interface 
profiles with the time sequences chosen in this review are at time t - 334, 392, 
431, 456, 487, 524, 554, 588, 600, 695, 784, and 800 s. It clearly shows that 


38 



there Is a series of asymmetric oscillations excited along the surface of 
sloshing dynamics governed liquid-vapor interface driven by asymmetric gravity 
gradient acceleration associated with slew motion. 

Comparison between Figures 9(A) and 10(A) shows the following results; (a) 
Components of gravity gradient acceleration in x and z directions are pointing 
toward the negative sides, which implies that the forces are pointing toward the 
southwestern direction. This indicates that the liquid is pushed toward the 
southwestern direction while the bubble is pushed toward the northeastern 
direction. These results are exactly shown in Figure 10(A). (b) Uneven and 
imbalance flow velocities toward the southwestern direction create similar uneven 
and imbalance pressure distribution reacting to bubble toward the northeastern 
direction. Because of extremely low surface tension coefficients of helium II 
between the liquid-vapor interface, a deformed irregular concave and convex- 
shaped oscillating bubble is created. 

Figure 11(A) shows time fluctuations of the locations of bubble mass 
centers of the fluids inside the dewar container due to sloshing dynamics driven 
by gravity gradient acceleration associated with slew motion. The values of 
bubble mass center fluctuations are (Ax c , Ay c , Az c ) — (11.8, 1.25, 18.4) cm. It 
shows Az c > Ax c > Ay c for bubble mass center fluctuations driven by gravity 
gradient acceleration associated with slew motion. 

The trend of the bubble mass center fluctuations, Az c > Ax c > Ay c , exactly 
reflects the values of the major driving forces of gravity gradient acceleration 
associated with slew motion in which it shows g r > g x > gy as that shown in 
Figure 9(A). 

Figures 12(A-a), 12(A-b), and 12(A-c) show the computed time variation of 
the fluctuations of slosh reaction forces acting on the dewar container driven 
by gravity gradient acceleration associated with slew motion in the y-axis with 
components along the x-, y-, and z-axes respectively. These figures show the 
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following results: (a) The values of slosh reaction force fluctuations are (AF X , 
AF y , AF t ) - (59.2, 4.11, 61.02) dynes; This clearly indicates AF Z > AF X >AF y . The 
maximum absolute values of slosh reaction force are Max (|F X |, |F y |, | F x | ) - 
(47.42, 4.11, 78.68) dynes. This also indicates |F X | > |F X | > |F y |. (b) The 
initial values of F x and F y start from zero, while that of F z starts from a 
nonzero value. This is due to the fact that the major driving force of gravity 
gradient acceleration associated with slew motion at the beginning of slew motion 
is equal to zero along the x- and y-axes, and is non-zero along the z-axis. (c) 
The characteristics of slosh reaction forces and their fluctuations are more 
likely to follow the trend of the major driving force of gravity gradient 
acceleration associated with slew motion in the y-axis as that shown in Fig. 11. 
Comparison between Fig. 12(A) and 9(A) reflects that the fluid system serves as 
a nonlinear modulator of spring mass and damper systems amplifying and damping 

the forces acting on the fluid flows and then responds back to the spacecraf t 58 ~ 
63 

Figures 12(B-a) , 12(B-b) , and 12(B-c) show the variations of slosh reaction 
moments acting on the dewar container due to bubble deformations driven by 
gravity gradient acceleration associated with slew motion with components along 
the x, y, and z axes, respectively. These figures show the following results: 
(a) The values of slosh reaction moment fluctuations are (AM^ AMy, AM Z ) - 
(770.6, 10784.5, 0.0004) dyne* cm. The maximum absolute values of slosh reaction 
moment are Max (IMJ, |My| , | M x | ) - (770.6, 8595.9, 0.0004) dyne- cm. (b) Because 
slew motion is in the y-axis, the magnitudes of both AMy and |My| are the 
maximum, (c) As slew motion is in the y— axis, the value of the moment arm along 
the y— axis is zero; and also, because the y— component of the major driving force 
of gravity gradient acceleration associated with slew motion in the y-axis is 
equal to zero, the magnitude of |F y | is near zero. Both |F y | ~ Ly - 0; thus, M* 
- L*F y - LyF x - 0. 


40 



Figures 12(C-a), 12(C-b) and 12(C-c) show time fluctuations of the moment 
arms of slosh reaction moments exerted on the dewar container due to bubble 
deformations driven by gravity gradient acceleration associated with slow motion 
in the y-axis for components along the x-, y- f and z-axes, respectively. These 
figures show the following results: (a) The values of moment arm fluctuation are 
(AL*, ALy, ALj) - (119.2, 11.7, 156.1) cm. The maximum absolute values of slosh 
reaction moment are Max ([LJ, |Ly|, 1 1^. | ) - (91.3, 11.7, 257.8) cm. This shows 
AL* > ALy > ALj and (L*) > |Ly| > |L,|. (b) Spacecraft slew motion s along the 
y-axis; thus, the value of Ly is the minimum, and that of My is the maximum, (c) 
Because L* is the function of FJiy, the value of My is the maximum, and that of 
F x is near the maximum. Thus, the values of L, are maximum. (d) The trend of 
the fluctuations of momentum arm is very much like that of the slosh reaction 
forces acting on the dewar container of the spacecraft. Also, the trend of the 
fluctuations of moment arm is very much like that of the major driving forces of 
gravity gradient acceleration associated with slew motion in the y-axis [see 
Figs. 12(C) and 9] . 

(VIII-B) Bubble Deformation Driven by G— Jitter Acceleration 
Associated With Slew Motion 

By using the mathematical formulations illustrated in Equations (2-3) to 
(2-11) and (2-18), one can numerically simulate spacecraft sloshing dynamics 
associated with spinning and/or slew motions depending upon the specific 
scientific missions assigned to the spacecraft. 

If slew motion operates with a range of 90° in 10 min (- 600 s) , the 
component of g-jitter acceleration along the (x, y, z) directions acted on the 
fluid mass located at (r, 6 , z) - (12 cm, n/ 2, 3 cm) is shown in Figure 9(B). 
This figure shows that the magnitude of slew motion induced acceleration is on 
the order of 10” 5 g^ 

The equilibrium shape of the liquid-vapor interface for a dewar with 70% 
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liquid-fill level under a residual gravity environment below 10“ 6 g 0 is a sphere. 
Figure 10(B) shows the time sequence evolution of the three-dimensional behavior 
of the interface oscillations driven by g-jitter acceleration associated with 
slew motion. It is shown in this figure that a time sequence evolution of 
liquid-vapor interface profiles at time t - 60.4, 258, 323, 354, 379, 403, 430, 
528, 572, 628, 663, and 800 s are illustrated. It clearly shows that there are 
a series of asymmetric oscillations excited along the surface of sloshing 
dynamics governed liquid-vapor interface driven by asymmetric g-jitter 
acceleration associated with slew motion 57-59 . 

For the convenience of explanation for the movement of bubbles at various 
locations in the figures of x-y, x-z and y-z planes, let us adopt the following 
conventional ways of direction illustration: (1) Choose positive-x, negative-x, 
positive-y, negative-y as the eastern, western, northern, and southern 
directions, respectively, in the x-y plane; (2) Choose positive-x, negative-x, 
positive— z, and negative-z as the eastern, western, northern and southern 
directions, respectively, in the x-z plane; and (3) Choose positive-y, negative— 
y, positive-z, and negative-z as the eastern, western, northern and southern 
directions, respectively, in the y-z plane. 

Figure 9(B) shows that g-jitter acceleration associated with slew motion 
is pointing toward the southeastern direction at the very beginning and then is 
pointing toward the southwestern direction at the end in the x-z plane, and also 
is pointing toward the southern direction in the y-z plane. These indicate that 
liquid is pushed toward the southeastern and then toward the southwestern 
directions (the bubble is pushed toward the northwestern and then toward the 
northeastern directions) while the liquid is pushed toward the southern direction 
(the bubble is pushed toward the northern direction) . 

Based on Figure 9(B), it shows that g-jitter acceleration associated with 
slew motion is pointing toward the eastern direction at the beginning and then 
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pointing toward the western direction at the end in the x-y plane. This implies 
that the liquid is pushed toward the eastern direction and then toward the 
western direction while bubble is pushed toward the western direction and then 
toward the eastern direction in the x-y plane. 

Comparison between Figures 9(B) and 10(B) illustrate some peculiar behavior 
of helium fluids with temperature below A-point (2.17 K) in which helium 
demonstrates a number of remarkable properties of superfluidity such as extremely 
low viscous and surface tension coefficients reacting to the disturbances driven 
by jitter acceleration associated with slew motion. It can be concluded as 
follows: (a) G-jitter acceleration associated with slew motion started at time 
t “ 0. However, bubble mass center fluctuations did not start to react to the 
driven force until t - 80 s for x c and t - 160 s for z c . In other words, bubble 
configuration was in perfect spherical-shaped at 80 s after g-jitter acceleration 
associated with slew motion was applied. (b) x— component g— jitter acceleration 
associated with slew motion was applied to the fluid element with positive value 
from t - 0 to 300 s; negative value from t - 300 to 600 s; and zero value after 
t — 600 s. However, westward movement of bubble continued from t — 80 to 450 s, 
and then switched to eastward movement t - 800 s which was 200 s after 
acceleration vanished. Obviously, there is a phase— lag between action of force 
and reaction of motion. (c) The z-component g-jitter acceleration associated 
with slew motion was applied to the fluid element with negative value from t — 
80 to 500 s and the zero value after t - 500 s. However, the northward movement 
of the bubble started at t - 160 s; bounced back from the wall at t - 530 s; and 
then the southward movement continued to t - 800 s which was 300 s after 
acceleration vanished. Obviously, the motion continued for a long period of time 
even after the applied force vanished due to extremely low viscosity of helium 
II fluids. (d) An intensive oscillation of bubble with a deformation of 
irregular concave and convex— shaped continued and sustained for several hundred 
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seconds after the applied force vanished due to extremely low surface tension 
coefficient for helium II fluids. 

Figure 11(B) shows time evolution of the locations of bubble mass centers 
of fluids inside the dewar container due to sloshing dynamics driven by g-jitter 
acceleration associated with slew motion. The values of bubble mass center 
fluctuations are (Ax c> Ay c , Az c ) - (28.9, 0.44, 30.2) cm. It shows Az c > Ax c > 
Ay c for bubble mass center fluctuations driven by jitter acceleration associated 
with slew motion. Behavior of bubble mass center fluctuations, shown in Figure 
11(B) are the exact reflection of the behavior of jitter accelerations, shown in 
Figure 9(B). 

Comparison between Figures 9(B), 10(B) and 11(B) illustrate some peculiar 
behavior of helium fluids with temperature below A-point (2.17 K) in which helium 
demonstrates a number of remarkable properties of superfluidity such as extremely 
low viscous and surface tension coefficients reacted to the disturbances driven 
by jitter acceleration associated with slew motion. It can be concluded as 
follows: (a) g-jitter acceleration associated with slew motion started at time 
t — 0. However, bubble mass center fluctuations did not start to react to the 
driven force until t - 80 s for x c and t - 160 s for z c . In other words, bubble 
configuration was in perfect spherical-shaped at 80 s after g-jitter acceleration 
associated with slew motion was applied, (b) x-component g— jitter acceleration 
associated with slew motion was applied to the fluid element with positive value 
from t - 0 to 300 s; negative value from t - 300 to 600 s; and zero value after 
t - 600 s. However, westward movement of bubble continued from t - 80 to 450 s, 
and then switched to eastward movement to t - 800 s which was 200 s after 
acceleration vanished. Obviously, there is a phase-shift between action of force 
and reaction of motion. (c) The z-component g-jitter acceleration associated 
with slew motion was applied to the fluid element with negative value from t - 
80 to 500 s and the zero value after t - 500 s. However, the northward movement 
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of the bubble started at t - 160 s; bounced back from the wall at t - 530 s; and 
then the southward movement continued to t - 800 s which was 300 s after 
acceleration vanished. Obviously, the motion continued for a long period of time 
even after the applied force vanished due to extremely low viscosity of helium 
II fluids. (d) An intensive oscillation of bubble with a deformation of 
irregular concave and convex-shaped continued and sustained for several hundred 
seconds after the applied force vanished due to extremely low surface tension 
coefficient for helium II fluids. 

Figures 13(A— a), 13(A— b) and 13(A— c) show the computed time variation of 
the fluctuations of slosh reaction forces exerted on the dewar container driven 
by g-jitter acceleration associated with slew motion in the y-axis with 
components along the x-, y— , and z-axes, respectively. These figures show the 
following results: (a) The values of slosh reaction forces are (AF X , AF y , AF t ) 
- (2.27, 0.04, 1.37) 10 3 dynes, it clearly indicates AF X > AF X > AF y ) . The 
maximum absolute values of slosh reaction forces are Max | F x | > |F X | > |F y | - 
(1.15, 0.02, 1.20) 10 3 dynes. It also indicates |F X | > |F X | > |F y |. (b) The 
initial values of F x , F y , and F x all start from zero value because there was no 
slew motion applied to the system at time « 0. (c) The magnitudes and 
fluctuations of F y are much smaller than those of F x and F z because the major 
driving moment of slew motion is along the y-axis. (d) The characteristics of 
slosh reaction forces and their fluctuations are ore likely to follow the trend 
of the major driving force of g-jitter acceleration associated with slew motion 
in the y-axis as that shown in Figure 9(B). Comparison between Figures 13(A) and 
9(B) reflects that the fluid systems serves as a damping modulator for the forces 
acting on the fluid flows and then responds back to the dewar container are due 
to the disturbances caused by the sloshing dynamics of the fluid system. 

Figures 13(B-a), 13(B-b), and 13(B-c) show the variations of slosh reaction 
moments exerted on the dewar container due to bubble deformations driven by g — 
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jitter acceleration associated with slew motion with components along the x, y, 
and z axes, respectively. These figures show the following results: (a) The 
values of slosh reaction moment fluctuations are (AM,, AMy, AMj) — (7.53, 422.79, 
0.0001) 10 3 dyne •cm. The maximum absolute values of slosh reaction moment are 
Max (IM*], I My- 1 • 1^*1) “ (4.31, 214.53, 0.0006) 10 3 dyne»cm. (b) As slew motion 
is in the y-axis , the magnitudes of both AMy and |My| are the maximum. (c) As 
slew motion is in the y-axis, the value of the moment arm along the y-axis is 
zero; and also as the y-component of the major driving force of g— jitter 
acceleration associated with slew motion in the y— axis is equal to zero, the 
magnitude of |F y | is near zero. Since both |F y | - Ly - 0, it leads M* - LyF, - 
L,F y - 0 and M, - L*Fy - LyF x - 0. 

Figures 13(C-a), 13(C-b) and 13(C-c) show the time fluctuations of the 
moment arms of slosh reaction moments exerted on the dewar container due to 
bubble deformations driven by g-jitter acceleration associated with slew motion 
in the y— axis for components along the x-, y— , and z— axes, respectively. These 
figures show the following results: (a) The values of moment arm fluctuations 
are (ALy, ALy, ALj) — (190.4, 56.1, 186.6) cm. The maximum absolute values of 
moment are ( 1 1^ | , |I*| , | Ly | ) - (95 .4 , 49 . 8 , 186 . 6) cm. It shows AL* > ALy > AL,) 
and |Lj| > |L*| > |Ly|. (b) As spacecraft slew motion is along the y-axis, the 
values of Ly are the minimum, and My are the maximum. (c) As Lj is the function 
of (FyM, - 0 because F y - 0) , the values of My are the maximum and F x is near 
the maximum. This makes the values of Ly the maximum. (d) The trend of the 
fluctuations of omentum arm is very much like that of the slosh reaction forces 
acting on the dewar container of the spacecraft. Again, the trend of the 
fluctuations of moment arm is very much like that of the major driving forces of 
g-Jitter acceleration associated with slew motion In the y-axis [see Figures 
13(C) and 9(B)]. 

IX. Rotating Bubble Subject Various Magnitudes of Gravity 
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Gradient and g-Jitter Accelerations 


The characteristics of gravity gradient and g-jitter accelerations, shown 
in Equations (2-14) and (2-18), respectively, are quite different. Gravity 
gradient acceleration shows the following characteristics: (a) acceleration acts 
on any fluid mass inside the container which increases two units of acceleration 
per unit of distance measured from the container mass center to the location of 
the fluid mass parallel along the radial axis from the container mass center to 
the center of the Earth; (b) the acceleration acts on the fluid mass which 
decreases one unit of acceleration per unit of the shortest distance measured 
from the location of the fluid mass to the radial axis along the vector from the 
container mass center to the center of the Earth 1 . For example, Figure 14(A) 
shows the time variation of gravity gradient accelerations for an assumed turn- 
around period of 1200 s with a container rotating speed of 0.1 rpm for components 
along (x, y, z) directions acted on the fluid mass located at (r, 0, z) — (40 cm, 
w/4, 10 cm). As the magnitude and direction of gravity gradient acceleration 
acting on each fluid mass is strongly dependent upon how far the location of the 
fluid mass deviates from the container mass center measured along the axis 
parallel to the vector r c which varies with respect to time, it shows that the 
gravity gradient acceleration acting on the fluid mass is different for fluid 
mass at different locations in the container. Figure 14(A) shows that the 
magnitude of gravity gradient acceleration is on the order of 1CT 7 g^ On the 
contrary to gravity gradient acceleration which exerts different acceleration on 
the fluid mass at different locations in the container, g— jitter acceleration 
drives the same acceleration on the fluid mass at different locations in the 
container. Figure 14(B) shows the time variation of g-jitter accelerations for 
a turn-around period of 1200 s with a container rotating speed of 0.1 rpm and a 
g-jitter frequency of 0.1 Hz for components along (x, y, z) directions acted on 
the fluid mass everywhere in the container. It is shown in Figure 9(A) that 
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gravity gradient accelerations are on the order of 10 -7 g Q for the spacecraft 
orbit chosen in this example, three ranges of background gravity, 10" 6 , and 10 -8 
go for g— jitter accelerations correspond to the magnitude of acceleration higher, 
and lower, respectively, than that of the gravity gradient acceleration acting 
on the fluid system of the spacecraft. 

The equilibrium shape of the liquid-vapor interface for a rotating dewar 
under a residual gravity environment below lO^go and rotating speed of 0.1 rpm 
is very much alike and is a doughnut configuration with a near circular kidney- 
shape cross-section based on the computation of the numerical algorithm developed 
in our earlier studies 7 " 9,20 . In this study, the combined gravity gradient and 
three ranges of g-jitter (10" 6 , and 10 -8 g 0 ) accelerations acting on the liquid- 
vapor interface oscillations have been reviewed and investigated 60-6 * . 

(IX— A) Combined Gravity Gradient and 10" 8 g 0 Background g-Jitter Accelerations 
Acting on Liquid-Vapor Interface Oscillations 

In this case, the combined effects of gravity gradient and g-jitter 
accelerations are completely dominated by the gravity gradient acceleration when 
the 10" 8 go and less background g-jitter accelerations are applied simultaneously 
with gravity gradient acceleration [which is on the order of 10 -7 go as shown in 
Figure 14(A)] for the excitation of slosh waves along liquid— vapor interface 
oscillations . 

Figure 15(A) shows the time sequence evolution of the liquid-vapor 
interface profiles driven by these combined accelerations. For the convenience 
of comparison, figures of liquid-vapor interface profiles with the same values 
of the time sequences chosen throughout this section are at time t - 191, 354, 
380, 431, 503, 603, 825, 980, 995, 1050, 1080 and 1200 s. It clearly shows that 
there is a series of asymmetric oscillations excited along the surface of liquid- 
vapor interface driven by asymmetric gravity gradient-dominated acceleration. 

Careful examination of the case of gravity gradient-dominated acceleration- 
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effected liquid-vapor interface oscillations, gravity gradient acceleration shown 
in Equation (2-12), indicate that there are greater negative components of 
acceleration longitudinal to the direction aligned with spacecraft mass center 
to the Earth center, and smaller positive components of acceleration transverse 
to this direction. As we indicated in Equation (2-2), angle varies with time. 
This phenomena shows that the gravity gradient acceleration exerted on the 
spacecraft is equivalent to the combination of time dependent force with turn- 
around direction and torsional moment acting on the spacecraft when it is 
orbiting around the Earth. It shows in Figure 15(A) that the deformation of the 
bubble is created by asymmetric torsional moment with a twisting force. 65 " 69 

Figure 15(A) shows the time sequence evolution of the three-dimensional 
dynamical behavior of the interface oscillations driven by gravity gradient- 
dominated acceleration. It indicates that the bubble (helium vapor) configura- 
tions change from axial symmetric to asymmetric profiles at a plane aligned with 
the vector of gravity gradient acceleration. The gravity gradient-dominated 
acceleration produces a combination of time dependent torsional movement with 
tidal motion of bubble oscillations in the rotating dewar when it is orbiting 
around the Earth. 

Figure 14(C) shows the time evolution of bubble mass center at (x c , y c , z c ) 
for three cases described above. The time history of both gravity gradient and 
g-jitter accelerations are plotted in Figure 14(B). Time evolution of bubble 
oscillations driven by gravity gradient-dominated, and g-j itter-dominated 
acceleration are illustrated in Figures 15(A), and (B) , respectively. 

Figure 14(C-a) shows time fluctuations of the locations of bubble mass 
centers of the rotating container due to slosh waves excited by gravity gradient- 
dominated acceleration. The values of bubble mass center fluctuations are (Ax c , 
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Ay c , Az c ) — (4.54, 2.83, 2.32) cm. It shows Ax c > Ay c > Az c for bubble mass center 
fluctuations driven by gravity gradient-dominated acceleration. 

As mentioned earlier, the combined effects of gravity gradient and g-jitter 
accelerations are completely dominated by gravity gradient acceleration when the 
10" 8 and less g-jitter accelerations are applied simultaneously with gravity 
gradient acceleration which drive slosh reaction forces and its moment 
fluctuations due to bubble deformations. 

Figure 16(A) shows the computed time variation of the fluctuations of slosh 
reaction forces exerted on the dewar container driven by gravity gradient- 
dominated acceleration. This figure shows the following results: (a) The values 
of slosh reaction force fluctuations are (AF X , AF y , AF X ) - (23.7, 18.3, 61.7) 
dynes, it clearly indicates AF X > AF y > AF X . The maximum absolute values of 
slosh reaction force are Max (|F X |, |F y |, | F x | ) - (17.4, 11.5, 32.4) dynes. It 
also indicates | F x | > | F x | > | F y | . (b) The initial values of F x and F y start from 
zero value while that if the F x starts from non— zero value. This is due to the 
fact that the shape of rotating bubble is symmetric with respect to the x and y 
axes which makes F x and F y equal to zero at time t - 0 while the resultant force 

the rotating axis of z is not equal to zero because of pre-existing gravity 
force acting along the z-axis at the very beginning, (c) The characteristics of 
slosh reaction forces and their fluctuations are very much following the trend 
of dominated gravity gradient acceleration , shown in Figure 7. In other words, 
gravity gradient acceleration as a major input force, modulated by sloshing 
dynamics of the fluid system in the dewar, activate on the dewar container with 
response force. In reality, the fluid system, modulated by sloshing dynamics, 
acts as a damper in the spacecraft to damp out the acting force of gravity 
gradient acceleration 43-47 . 

Figures 16(B-a), 16(B-b), and 16(B-c) show variations of slosh reaction 
torque exerted on the dewar container due to slosh waves driven by gravity 
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gradient-dominated acceleration along the x, y, and z axes, respectively. The 
values of slosh reaction torque fluctuations are (AMy, AMy, AM*) - (761.9, 667.3, 
0.06) dyne* cm. The maximum absolute values of slosh reaction torques are Max 
( |My| , |My|, |M*|) - (392 . 9 , 343 .4, 0. 50) dyne* cm. It shows AMy > AMy > AM* and 
|My| > | My | > |M*|. As Mjj, My, M* are determined by the Factors (F x , L*, F XI Ly) , 
(F x , Ly, F y , L*) and (F y , Ly, F x , Ly) , respectively, it was illustrated in Figures 
32 and 34 that the magnitudes of force fluctuations give |F X | > |F X | > | F y | and 
that of moment arm fluctuations give IL*) > |Ly| > |Ly| in which time averages 
of Ly and Ly were approaching zero. This explains why the moment of M* - LyF y - 
LyF x - 0 even though there is a rolling moment with angular velocity of 0.1 rpm 
in which the induced rolling moment exerted by the viscous rotating fluids is 
small in comparison with the induced pitch and yaw moments due to asymmetric 
slosh wave excitation. 

Figure 16(C) shows time fluctuations of moment arms of slosh reaction 
torques exerted on the dewar container due to bubble deformations excited by 
gravity gradient-dominated acceleration. The values of moment arm fluctuations 
are (ALy, ALy, AL*) - (22.8, 19.0, 45.5) cm. The maximum absolute values of 
moment arm are Max (|LJ # |Ly|, | L* | ) - (12.2, 10.2, 103.4) cm. It shows ALy > 
ALy > AL* and |Ly| > |Ly| > | L* J . This trend is very much similar to the driving 
force of gravity gradient acceleration, shown in Figure 14(A) . 

(IX-B) Combined Gravity Gradient and 10~ 6 g<, Background g-Jitter 
Accelerations Acting on Liquid— Vapor Interface Oscillations 

In this case, the combined effects of gravity gradient and g-jitter 
accelerations are completely dominated by the g-jitter acceleration. 

Figure 15(B) shows the time sequence evolution of the dynamical behavior 
of the liquid— vapor interface oscillations driven by g-j itter— dominated 
acceleration of 0.1 Hz low g-jitter frequency. It clearly shows that there Is 
a series of asymmetric oscillations excited along the surface of liquid-vapor 
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interface driven by asymmetric g-j i tter-dominated acceleration. 

Careful examination of the case of g-jitter acceleration-dominate liquid- 
vapor interface oscillations indicates that there is a sinusoidal oscillation 
longitudinal to the direction aligned with spacecraft mass center to Earth 
center. As indicated in Equation (2-2), the angle varies with time. This 
phenomena shows that the g-jitter acceleration exerted on the spacecraft is 
equivalent to time-dependent oscillatory forces which push the bubble in the 
combined directions of down-and-up [see z-component of g— jitter acceleration 
shown in Figure 14(B)] and leftward-and-rightward [see x-and-y components of g- 
jitter accelerations shown in Figure 14(B)] as the bubble is rotating with 
respect to the spacecraft rotation axis. The down-and-up motion of the bubble 
is due to the fact that z-coraponent g-jitter acceleration decreases from the 
maximum downward direction to the minimum downward direction, during the time 
period between t - 0 to t - 300 s, which pushes the liquid downward and the 
bubble upward; while the g-jitter acceleration changes to its downward direction 
and magnitude from the minimum upward to the maximum upward, and then back to the 
minimum upward, during the time period between t - 300 s to t - 900 s, which 
pushes the liquid upward and the bubble downward; in the meanwhile, the g— jitter 
acceleration changes to its downward direction and magnitude during the period 
from t - 900 s to 1200 s which repeats the half cycle story of the time period 
from t - 0 to t - 300 s. The leftward-and-rightward oscillations of bubble are 
the results of positive and negative directions of g-jitter accelerations in the 
x- and y-components as that illustrated in Figure 14(B). The positive direction 
g-jitter acceleration implies that the liquid is pushed rightward and the bubble 
is driven leftward. The time variations of g-jitter accelerations, shown in 
Figure 14(B) , characterize the leftward and rightward oscillations of the bubbles 
shown in Figure 15(B). 

In addition to the 0.1 Hz low frequency j itter— dominated acceleration, 1.0 
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Hz medium frequency and 10 Hz high frequency jitter acceleration driven liquid- 
vapor interface oscillations have also been investigated. The characteristics 
of these three ranges of g-jitter frequencies have been well-documented in the 
earlier studies 10 " 20 . Because of page limitation, oscillations driven by the time- 
dependent variations of g-j it ter— dominated acceleration with medium and high 
frequencies will not be shown in this review. It is also worthwhile to mention 
that the natural frequency of a cryogenic helium container is much closer to the 
low frequency g-jitter acceleration than that of the medium and high frequencies. 
Results show that lower frequency g-j itter-dominated acceleration contribute more 
in the driving of asymmetric profiles of the interface than the higher frequency 
g-j itter-dominated acceleration. These results are agreeable with the earlier 
studies 10 " 20 . 

Figure 14(C-b) shows time fluctuations of the locations of bubble mass 
centers of the rotating container due to slosh waves excited by asymmetric 
gravity jitter-dominated acceleration. The values of bubble mass center 
fluctuations are (Ax c , Ay c , Az c ) - (17.28, 12.2, 33.83) cm. A comparison of 
Figures 14 (C— a) and 14(C-b) for bubble mass center fluctuations driven by gravity 
gradient-dominated, and g-j itter-dominated accelerations, provides the following 
conclusions: (a) Torsional moment— and twisting force-equivalent gravity 
gradient— dominated acceleration exerted on the rotating dewar container produces 
smaller values of magnitudes and fluctuations of angular momentum and fluid 
moment than that driven by g-j itter-dominated acceleration. (b) Dynamics of 
bubble (liquid— vapor interface) driven by torsional moment— and— twisting force- 
equivalent gravity gradient acceleration produces bubble mass center fluctuations 
of Ax c > Ay c > Az c while up and down oscillations of bubble driven by g-j itter- 
dominated acceleration induces bubble center fluctuations of Az c > Ax c > Ay c . 
(c) A comparison of the values of bubble mass center fluctuations show that 
fluctuations of bubble mass center induced by g-j itter-dominated acceleration is 
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greater than that induced by gravity gradient-dominated acceleration. (d) 
Fluctuations of both x c and y c start from zero while that of the z c starts from 
a non-zero value at the middle point of the height of partially filled liquid 
container. (e) Fluctuation magnitudes of x c , y c and z c associated with higher 
frequency jitter are lower than those of the frequencies associated with lower 
jitter frequency for bubble mass center fluctuations driven by g-j itter-dominated 
acceleration. 

Figure 17(A) shows the computed time variation of the fluctuations of slosh 
reaction forces exerted on the dewar container driven by g—j itter-dominated 
acceleration with frequency of 0.1 Hz. The values of slosh reaction force 
fluctuations are (AF X , AF y , AF Z ) — (383.0, 315.5, 619.7) dynes and the maximum 
values of slosh reaction force are Max ( | F x | , |F y |, | F z | ) - (192 . 2 , 172.3, 359.5) 
dynes. It shows AF X > AF y > AF Z and |F X | > |F y | > |F Z |. Comparison of input and 
response of forces between Figures 14(B) and 17(A) show that these two sets of 
curves are very much similar in characteristics and trends with smaller 
amplitudes and fluctuations shown in the response curve. In other words, it 
indicates that the fluid system, with modulation of sloshing dynamics, acts as 
a damping machine in input and response circuit. Comparison between Figures 
16(A) and 17(A) for slosh reaction force fluctuations driven by gravity gradient- 
dominated acceleration and g-j itter-dominated acceleration show that gravity 
jitter-dominated acceleration are responsible for the excitation of greater 
values and magnitudes than that driven by gravity gradient-dominated accelera- 
tion. 

Figures 17(B-a), 17(B-b), and 17(B-c) show variations of slosh reaction 
torque exerted on the container wall due to bubble deformations excited by g- 
j itter-dominated acceleration with frequency of 0.1 H along t ex, y, and z axes, 
respectively. The values of slosh reaction torque fluctuations are (AM X , AMy, 
AM Z ) — (1612.2, 1957.9, 0.10) dyne* cm and the maximum absolute values of slosh 
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reaction torque are Max (IM*), |My|, |M X |) - (1013.7, 1139.3, 0.52) dyne*cm. It 
shows AMy > AM* > AM X and |My| > | | > |M X |. Characteristics of the fluctuations 
of slosh reaction torques driven by g-j Itter-dominated acceleration are similar 
to that drawn for the fluctuations of slosh reaction torque driven by g-jitter- 
dominated acceleration. However, g-j itter-dominated acceleration excited values 
and fluctuations of slosh reaction torque exerted on the dewar container have a 
greater value and magnitude than that driven by gravity gradient-dominated 
acceleration. 

Figure 17(C) shows time fluctuations of moment arm of slosh reaction torque 
exerted on the dewar container due to bubble deformations excited by asymmetric 
g—j itter-dominated acceleration. The values of moment arm fluctuations are (AL*, 
ALy, AI^) — (4.6, 3.9, 5.3) cm. The maximum absolute values of moment arm are 
Max (|L*|, |Ly|, 1 4 1 ) - (2.4, 3.3, 77.8) cm. It shows AL* > ALy > AL* and | 4 | 
> 141 > 141 • Comparison of gravity gradient-dominated and g-jitter dominated 
accelerations driven fluctuations provide the following conclusions: (a) 
Torsional moment— and twisting force— equivalent gravity gradient— dominated 
acceleration exerted on the rotating dewar container produces smaller maximum 
values and fluctuations of slosh reaction forces and torques than that driven by 
g—j itter-dominated acceleration, (b) Dynamics of bubble (liquid— vapor interface) 
driven by torsional moment-and twisting force-equivalent gravity gradient- 
dominated acceleration produces greater fluctuations and magnitude of slosh 
reaction moment arm than that produced by up and down oscillations of bubble 
driven by g-j itter-dominated acceleration. (c) Input and response relation 
between combined gravity gradient/g-j itter accelerations and slosh reaction 
forces acting on the dewar container through the modulation of sloshing dynamics 
show that the characteristics and trends are similar. However, the fluid system 
acts a damper which, in a realistical sense, helps damping out the spacecraft 
acceleration activated by the orbital forces* 2 ”* 5 . 
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X Discussion and Conclusions 


Mathematical formulations of the possible accelerations which may exert on 
the fluid systems of spacecraft in microgravity have been reviewed. Dynamics of 
non-inertial frame fluid systems governing equations subject to initial and 
boundary conditions applicable to bubble dynamics in microgravity have been 
discussed. Bubble deformations and oscillations driven by (a) Spin-up from rest 
with and without a completion of wrapping around the dewar well, (b) Spin— down 
from complete wrapping of rotating dewar with steady state, (cc) rotating bubble 
subject to various magnitudes of gravity gradient and g— jitter accelerations 
associated with slew motion, (e) bubble deformations driven by orbital 
accelerations in rotating dewar with and without baffle, and (f) rotating bubble 
subject to lateral and axial impulses , have been reviewed and discussed. 

Bubble mass center fluctuations due to bubble deformations and oscillations 
driven by various magnitudes of orbital accelerations have caused major problems 
in spacecraft orbital and attitude controls 34 * 47 . This problem has been discussed 
precisely in this review. 

Slosh reaction forces and torques fluctuations acting on the dewar due to 
bubble deformations and oscillations driven by the orbital accelerations have 
been another source of problems induced in spacecraft controls 34 * 47 . Precise 
mathematical formulations of the computations of these slosh reaction and torques 
have been discussed. Accurate computation and prediction of these slosh reaction 
forces and torques are essential for the development of orbital and attitude 
control techniques of spacecraft 34 * 47 . 

Similarity rules have been used widely in the prediction of bubble 
oscillations in dewar container with different geometry and various physical 
parameters. With help of Hung et al 37 , we can extent the discussion covered in 
this review to various cases with different geometry. 

For the purpose to truly validate computation code developed, it is 
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urgently needed to carry out sloshing experiment in microgravity, in time 
evolution of both bubble deformation and slosh reaction forces and torque. 
Computational programs developed for slosh dynamics so far is quite complicated 
and is not user friendly. It is suggest to spend more effort in developing 
numerically efficient computational code to be used for future spacecraft real 
time operation. 
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Figure Captions 


Figure 1 


Figure 2 
Figure 3 


Figure 4 
Figure 5 


Figure 6 
Figure 7 
Figure 8 


AXAF-S spacecraft coordinate systems with azimuth angle from 
spacecraft mass center to the center of the Earth. Coordinate (x" , 
y" , z") for slew motion and coordinate (x, y, z) for fluid mechanics 
computation. 

Computational algorithm for sloshing dynamics. 

(A) Time sequence evolution of three-dimensional cryogenic helium II 
bubble fluctuations for spin-up rotating dewar. Liquid filled 
level-95%, w - 0.4 rpm, and We - 226.62. (B) Various equilibrium 
profiles of horizontal cross section at height z-72.5 cm. 

Time sequence evolution of cryogenic helium II bubble mass center 
fluctuations for spin-up rotating dewar with and without completion 
of bubble wrapping around dewar well. 

Final equilibrium bubble profiles in r-z plane and 6 - 0° and 180°, in 
r-0 plane at height z - 72.5 cm and in three-dimensional with liquid 
level of 85% (A) w - 0.15 rpm, We - 31.87, t - 480 s, (B) w - 0.05 rpm, 
We - 3.54, t - 1600 s. 

(A) Three-dimensional configuration of container with baffle, (B) 
distribution of grid points in the radial-axial plane, and (C) in the 
radial-circumferential plane. 

Time evolution of three-dimensional bubble oscillations for time 
between 225 and 427 s for rotating dewar in response to lateral impulse 
without and with baffle. 

(A) Time sequence evolution of fluid mass center fluctuations in 
response to lateral impulse for container with and without baffle, 
(a) Fluid mass center fluctuations in x-axis, (b) Fluid mass center 
fluctuations in y— axis , and (c) Fluid mass center fluctuations in 
absolute values of radial direction. (B) Time sequence evolution of 
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fluid mass center fluctuations in response to axial impulse for 
container with and without baffle. 

Figure 9 (A) Time variation of AXAF-S spacecraft gravity gradient acceleration 

acting on fluid mass located at (r, 0, z) - (12cm, *r/2 , 3 cm) for 90° 
slew motion in 600 s along the y"-axis and orbital period of 97.6 min. 
(B) Time variation of AXAF-S spacecraft jitter acceleration associated 
with slew motion acting on fluid mass located at (r, 0, z) — (12 cm, 
n/2 , 3 cm) for 90° slew motion in 10 min. along the y"-axis and orbital 
period of 97.6 min. 

Figure 10 (A) Time sequence evolution of three-dimensional bubble oscillations 
for dewar driven by gravity gradient acceleration associated with slew 
motion in the y"-axis. (B) Time sequence evolution of three-dimensional 
bubble oscillations for dewar driven by g-jitter acceleration 
associated with slew motion in the y"-axis. 90° slew motion in 600 s 
is applied to the spacecraft operation. 

Figure 11 (A) Time sequence of the AXAF-S spacecraft fluid moment fluctuations 
due to sloshing dynamics driven by gravity gradient acceleration 
associated with slew motion in the y-axis. (B) Time sequences of the 
AXAF— S spacecraft bubble mass center fluctuations due to the sloshing 
dynamics driven by g-jitter acceleration associated with slew motion 
in the y-axis. 90° slew motion in 600 s is applied to spacecraft 
operation. 

Figure 12 Time sequence of the fluctuations of (A) slosh reaction forces (B) 
slosh reaction moment, and (C) slosh reaction moment arm exerted on 
the dewar due to bubble deformations driven by gravity gradient 
acceleration associated with slew motion. 

Figure 13 Time sequence of the fluctuations of (A) slosh reaction forces (B) 
slosh reaction moment, and (C) slosh reaction moment arm exerted on 
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the AXAF-S dewar due to sloshing dynamics driven by g-jitter accelerat- 
ion associated with slew motion. 

Figure 14 (A) Time variation of gravity gradient acceleration acting on fluid 
mass located at (r,0,z)-(4O cm, */4, 10 cm) for turn-around period of 
1200 s with rotating speed of 0.1 rpm. (B) Time variation of jitter 
acceleration acting on the fluid system under background gravity of 
10~ 8 , 10~ 7 and 10“ 6 g ol rotating speed of 0.1 rpm, turn-around period 
of 1200 s and jitter frequency of 0.1 Hz. (C) Time sequence of bubble 
mass center fluctuations caused by slosh wave excitation under gravity 
turn-around time of 1200 s rotation speed of 0.1 rpm. 

Figure 15 Time sequence evolution of three-dimensional liquid-vapor interface 
oscillations for rotating dewar driven by (A) combined gravity gradient 
and jitter accelerations with magnitude of 10“ 8 go and (B) combined 
gravity gradient and jitter accelerations of 0.1 Hz frequency, and 
with magnitude of 10“ 6 g Q . 

Figure 16 Time sequence of the fluctuations of (A) slosh reaction forces (B) 
slosh reaction moment and (C) slosh reaction moment arm exerted on 
the dewar container due to the bubble deformations driven by combined 
gravity gradient and g-jitter accelerations with 10 -8 g 0 under 
rotating speed of 0.1 rpm and gravity turn-around time of 1200 s. 

Figure 17 Time sequence of the fluctuations of (A) slosh reaction forces (B) 
slosh reaction moment, and (C) slosh reaction moment arm exerted on 
the dewar due to the bubble deformations driven by combined gravity 
gradient and g-jitter accelerations with 10” 7 go, rotating speed of 0.1 
rpm, gravity turn-around time of 1200 s and jitter frequency of 0.1 Hz. 
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(A) Gravity Gradient Acceleration ( B ) Gravity Jitter Acceleration 
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THREE DIMENSIONAL LIQUID AND VAPOR INTERFACE 
The Effect of Gravity Jitter and Gravity Gradient 


















Slosh Reaction Force, Torque. and Moment' Arm 
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